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ABSTRACT 

The postcollapse structure of objects which form by gravitational condensation out 
of the expanding cosmological background universe is a key element in the theory of 
galaxy formation. Towards this end, we have reconsidered the outcome of the nonlin¬ 
ear growth of a uniform, spherical density perturbation in an unperturbed background 
universe - the cosmological “top-hat” problem. We adopt the usual assumption that 
the collapse to infinite density at a finite time predicted by the top-hat solution is in¬ 
terrupted by a rapid virialization caused by the growth of small-scale inhomogeneities 
in the initial perturbation. We replace the standard description of the postcollapse 
object as a uniform sphere in virial equilibrium by a more self-consistent one as a trun¬ 
cated, nonsingular, isothermal sphere in virial and hydrostatic equilibrium, including 
for the first time a proper treatment of the finite-pressure boundary condition on the 
sphere. The results differ significantly from both the uniform sphere and the singular 
isothermal sphere approximations for the postcollapse objects. The virial temperature 
which results is more than twice the previously used “standard value” of the postcol¬ 
lapse uniform sphere approximation but 1.4 times smaller than that of the singular, 
truncated isothermal sphere approximation. The truncation radius is 0.554 times the 
radius of the top-hat at maximum expansion, and the ratio of the truncation radius 
to the core radius is 29.4, yielding a central density which is 514 times greater than 
at the surface and 1.8 x 10 4 times greater than that of the unperturbed background 
density at the epoch of infinite collapse predicted by the top-hat solution. For the 
top-hat fractional overdensity Sl predicted by extrapolating the linear solution into 
the nonlinear regime, the standard top-hat model assumes that virialization is instan¬ 
taneous at Sl = S c = 1.686, i.e. the epoch at which the nonlinear top-hat reaches 
infinite density. The surface of the collapsing sphere meets that of the postcollapse 
equilibrium sphere slightly earlier, however, when Sl = 1.52. These results will have 
a significant effect on a wide range of applications of the Press-Scliechter and other 
semi-analytical models to cosmology. 

We discuss the density profiles obtained here in relation to the density profiles 
for a range of cosmic structures, from dwarf galaxies to galaxy clusters, indicated by 
observation and by N-body simulation of cosmological structure formation, includ¬ 
ing the recent suggestion of a universal density profile for halos in the Cold Dark 
Matter (CDM) model. The nonsingular isothermal sphere solution presented here pre¬ 
dicts the virial temperature and integrated mass distribution of the X-ray clusters 
formed in the CDM model as found by detailed, 3D, numerical gas and N-body dy¬ 
namical simulations remarkably well. This solution allows us to derive analytically 
the numerically-calibrated mass-temperature and radius-temperature scaling laws for 
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X-ray clusters which were derived empirically by Evrard, Metzler and Navarro from 
simulation results for the CDM model. 

Key words: cosmology: theory - dark matter - galaxies: clusters: general - galaxies: 
formation - galaxies: haloes - galaxies: kinematics and dynamics 


1 INTRODUCTION 


The theory of structure formation in cosmology must answer the question of how gravitational instability amplified initially 
small-amplitude density fluctuations in an expanding background universe until the nonlinear structure we see today, from 
galaxies to clusters of galaxies, emerged. It must explain the number and frequency of these objects, their spatial clustering and 
relative motions, their internal structure, and the evolution of these properties. The current view is that this structure arose 
by the complex, hierarchical build-up which resulted from a primordial spectrum of Gaussian-random noise which reached 
nonlinear amplitude some time after matter began to dominate the total energy density of the universe. The exact nature of 
this structure formation is not amenable to an analytical solution. As a result, much effort has gone into large-scale numerical 
simulation by N-body techniques and, more recently, the coupling of N-body and hydrodynamics methods. For many purposes, 
however, analytical and semi-analytical approximations are possible and serve both as a guide to understanding the complex 
numerical results and as a crucial means of extrapolating beyond them in order to interpret the predictions of cosmological 
models in comparison with the observed universe. 

The nonlinear evolution of cosmic density fluctuations is often approximated by a model in which the initial linear 
perturbation is an isolated, uniform sphere outside of which the matter is unperturbed—the “top-hat” perturbation — for 
which an exact solution is possible (e.g. Peebles 1980, Peebles 1993, Padmanablian 1993; and references therein). For many 
purposes, a critical part of the application of this top-hat model involves the assumption that the collapse to infinite density 
predicted by the solution is interrupted by a rapid relaxation to virial equilibrium at a finite density, as a result of the 
growth of initia lly small inhomogeneities in the density distribution. For example, many applications of the Press-Schechter 
approximation (Press and Schechter 1974) for the rate of formation of objects of a given mass per unit time due to the growth 
of Gaussian-random-noise density perturbations involve this assumption of postcollapse virialization (see, for example, Monaco 
1997 for a review and references). 

We have reconsidered the postcollapse virial equilibrium of the top-hat model. We improve upon the standard assumption 
of a uniform sphere for the virialized final state by a self-consistent model involving a nonsingular, truncated isothermal sphere 
(“Tig”) j n hydrostatic and virial equilibrium which takes proper account of the finite-pressure boundary condition on the 
sphere. The results differ significantly from those of the standard uniform sphere (“SUS”) approximation for the final virialized 
object. They also differ significantly from those of the approximation of a singular isothermal sphere (“SIS”) for the final 
object. 

An outline of the paper is as follows. In § b we describe the top-hat model and the standard approach to the final 
virialized object in the SUS approximation. In § tj we describe the hydrostatic equilibrium solutions for isothermal spheres. 
An application of the virial theorem to characterize the equilibrium of a truncated isothermal sphere, including a proper 
account of the finite boundary pressure, is described in § |S|. A minimum-energy argument which selects a unique TIS solution 
for a given boundary pressure is presented in § g This solution is used to specify the postcollapse TIS which results from 
any given top-hat density perturbation in § g. An independent calculation of this postcollapse TIS solution is presented in 
§ g based on the self-similar infall solution of Bertschinger (1985), with results which agree closely with the minimum-energy 
argument, thereby confirming it. Our results are summarized and discussed in § H. Our summary in § 8.1 includes a complete 


prescription, with scaling laws and an approximate fitting formula, for using our minimum-energy TIS solution in combination 
with the top-hat perturbation solution to describe the cosmological formation and postcollapse equilibrium structure of any 
object as a function of its mass and collapse epoch. In § |S.2| , we compare the predictions of this TIS solution with those of the 
SUS and SIS approximations. A comparison of the TIS density profile with existing results of numerical simulations of cosmic 
structure formation and with observations for a range of cosmic structures from dwarf galaxies to galaxy clusters is made in 
§ This includes a discussion of the minimum-energy TIS solution in relation to the suggestion by Navarro, Frenk, and 
White (1996, 1997) of a universal form for the density profiles of haloes of collisionless dark matter formed in a hierarchical 
c luste ring model like the Cold Dark Matter (CDM) model, based upon the results of numerical N-body simulation. Finally, in 
§ 8.4, we compare our predicted virial temperature with the results of numerical gas dynamical simulations of X-ray cluster 
formation in the CDM model, and derive analytically the mass-radius-temperature scaling laws for X-ray clusters determined 
empirically by Evrard, Metzler, and Navarro (1996) from gas dynamical and N-body simulations of the CDM model. 
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2 THE STANDARD APPROACH TO THE TOP-HAT MODEL 

All the results in this section are well known, but due to their importance in what follows we shall summarize them here. 


2.1 Before collapse: the exact nonlinear solution 


The spherical top-hat model, a spherical perturbation of uniform overdensity ( Gunn fe Gott 1972 ), affords considerable insight 
into the dynamics of gravitational growth of cosmic structure, while still having an exact, analytical or semi-analytical, solution 
in certain cosmological models (cf. Peebles 1980). In what follows we shall consider top-hat perturbations only in the Einstein- 
de Sitter model. 

Based on its high degree of symmetry, the collapse of a top-hat density perturbation can be described in detail solely in 
terms of its overdensity with respect to the background. From the Birkhoff theorem, the equation describing this evolution 
is identical to the Friedmann equation for a universe with the corresponding mean density. The linearized version of this 
equation describes the linear evolution of the density perturbation 6z,(t) and has a growing solution Sl = Si(t/ti) 2 ^ 3 . Using 
this solution we can write the equation for the density contrast in terms of 5n(t), whose solution is independent of the 
initial density contrast Si and the corresponding time ti and hence describes the general, scale-invariant collapse solution. In 
parametric form the solution reads (e.g. Padmanabhan 1993) 


9 (9 — sin9) 2 
d ^2(l-cos0) 3 


&-1 (ir <*->»<’■ 


(i) 


The critical density contrast S c is defined as the linear solution extrapolated to the epoch at which the nonlinear solution 
predicts an infinite density. As described above, this effectively is the time of collapse in scale-free units. The values obtained 
at turnaround (9 = n, indicated by subscript “m” to denote “maximum expansion”) and at collapse time (9 = 2n, denoted 
by subscript c) are: 


3 /37t\ 2 / 3 Q7t 2 3 /37t\ 2 / 3 

SLm = - J « 1.0624, 5m = ——1 «4.5517, 5 C = - y —J « 1.6865, S(S C ) = oo. (2) 

Henceforth, we shall refer to z co i\ as the redshift which corresponds to the epoch of infinite collapse, at which 5 = oo, at time 

tcoll- 


2.2 After collapse: uniform sphere in virial equilibrium 


As shown above, a perfectly symmetric top-hat collapse results in a singularity at 5l = S c , which does not lead to the formation 
of virialized structures. To overcome this problem, it is usually assumed that the actual collapse is slightly inhomogeneous, 
and therefore the top-hat does not collapse to infinite density, but instead, by means of processes like violent relaxation, 
relaxes to form a static, virialized structure. Assuming that the total energy is conserved during the collapse, we can connect 
the initial top-hat to the final state as follows. 

The standard approach is to assume that the collapse of the top-hat to infinite density is interrupted by a rapid equi¬ 
libration at Sl = Sc which results in another uniform sphere in virial equilibrium. We solve for the final radius r v ir of the 
virialized sphere by applying the virial theorem to the final state and expressing the conserved total energy E of the sphere 
in terms of the radius r m at maximum expansion. At the point of maximum expansion (which always exists for fl > 1) the 
sphere is cold and at rest, so its energy is entirely gravitational potential energy which, for uniform sphere of mass Mo and 
radius r m is just 


E = Wm = 


3 GMp 
6 r m 


(3) 


After the system settles down to a virial equilibrium its potential energy is 


ITvir = - 


3 GA/ 0 2 
5 Tvir 


(4) 


According to the virial theorem 


0 = 2 K + W, 


(5) 


where K = Uth + 2km, Uth and Tki n are the thermal and kinetic energy, respectively. Together 
E = K + W, this implies E = —W/2, and, therefore, the total energy is related to the virial 

E 3 GMl 

10 r v ir 


with the conservation of energy, 
radius according to 


( 6 ) 


Equating the total energy E in equations (pi) and (pi) yields 
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( 7 ) 


This is the well-known result for the size of the virialized object in the SUS approximation. The virial temperature and velocity 
are derived as follows. The kinetic energy in the virialized state is 


K v ir = E- Wvir = 


3 GMl 


10 TVir 

This kinetic energy is the energy of internal motions only (Tkin =0). To obtain the equivalent temperature we express the 
thermal energy 

3 k B T . . 

U t h — r- M 0 , (9) 

2 m 

where m is the mean mass per gas particle. If m B is the mass of a hydrogen atom, then m = pm B where p is the mean 
molecular weight. Therefore the equivalent virial temperature is 
T _ 1 GM 0 m _ 2 GM 0 m 
5 1’ j3 T’ vir 5 k H 

Henceforth, we shall refer to the virial temperature in equation © as Tsus- For the case of a collisionless gas, we replace 
the virial temperature above by the virial velocity dispersion, 

■■ 2x k B T 


( v ~. 

3 


( 11 ) 


According to this approach, the properties of the final virialized uniform sphere are related to those of the top-hat 
perturbation which created it as follows, as summarized by Padmanabhan( 1993| ) , except we use h = Ro/(100kms _1 Mpc -1 ) 
rather than ho.s = 2 h, and we include the explicit dependence on the mean molecular weight p, where p = 0.59 (1.22) for an 
ionized (neutral) gas of H and He with [He]/[H] = 0.08 by number: 

( M V /3 Ci. — lfi — ^/3 


r m = 338.3 


V 1O 12 M 0 


(1 + Zcoii) h 7 kpc, 


{r v ir)sus = 169.2 ^ 10 i 2 j\/ e ^ ( : + 2co11 ) 2/3k P c > 

/ \ 2/3 

Tsus = 6.160 x 10 5 p f 1Q1 J (1 + 0coii)A 2 / 3 K, 


( 12 ) 

(13) 

(14) 


{v c )sus = 5 1/2 a = 



159.4 


M \ 
10 12 Mq) 


1/3 

(1 + a C oii) 1/2 h 1 ' /3 kms _1 , 


(15) 


(. p)sus = 187T 2 (1 + 0coll) 3 pbO, 


(16) 


where v c = ( GM/r v i T ) J ^ 2 is the circular velocity at r v ir- 

This approach is commonly used to calculate the virial temperature and radius. However it is not a realistic model, 
since it assumes isothermality, uniform density, and the absence of external pressure, which are clearly incompatible. A more 
realistic approach is to assume a final state of hydrostatic equilibrium (or the corresponding state with an isotropic Maxwellian 
velocity distribution in the case of collisionless particles), which we will do in the next section. 


3 ISOTHERMAL SPHERES 


The final virialized object is decoupled from the expanding cosmological background from which it condensed. Hence, when 
we describe it as an isothermal sphere in hydrostatic equilibrium, we do so in the usual non-cosmological way (cf. Binney and 
Tremaine 1987). The hydrostatic equilibrium equation, Vp = pg, in the case of spherical symmetry becomes 


k B dp GM(r) 

T —— = -pg = -p - - — 

m dr r z 


(17) 


where M(r) is the mass inside radius r. Multiplying equation (|I^) by r 2 m/pk B T, and taking the derivative with respect to 
r, we obtain 


d 

dr 


; ri(lnp) 

dr 


Gm 2 

= — 47t-—— r p. 
k B T H 


(18) 


Let us consider the case of collisionless particles, too. The Poisson equation for the gravitational potential in the case of 
spherical symmetry is 
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Id/ 2 d'& 
r 2 dr V dr 


)- 


47t Gp. 


The equilibrium velocity distribution of the particles is a Maxwellian distribution given by distribution is 


f(v) = 


Po 


(27rcr 2 ) 3 / 2 


exp 


VE< -v 2 /2 


(19) 


( 20 ) 


where po is the central density if we take \&(r = 0) = 0, and a is the one dimensional velocity dispersion. After integrating 
over velocity we obtain 


P = J f(v)dv = poe !/a 

Therefore 


T = rr In 2- 
\Po 

which we substitute in equation (|lo|), obtaining 
2 d(lnp) 


d 

dr 


dr 


4tt 2 
=- 2 Gpr . 


By calculating the mean squared velocity we obtain 
(' V 2 ) = 3cr 2 . 

The equivalent temperature can be calculated from 

(v 2 ) = 3 kgT 
2 2 m 

obtaining for cr. 

2 ksT 

a = -. 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


A comparison of equation (j2^) with equation @ using equation (|2^) shows they are identical. Hence, the structure of an 
isothermal self-gravitating fluid sphere is identical with the structure of a collisionless 3D spherical system in equilibrium. 

To make equation (E3|) nondimensional, we introduce new variables 


- _ P , _ r 

Pi C 5 

po r o 

where po is the central density, and we choose 
ro = cr /(AnGpo) 1 ^ 2 . 

In these variables equation (p3|) becomes 

A 2d(lnp)' 


d_ 

d( 


C 


d( 


= -PC- 


We must solve equation ([29|) with boundary conditions 
dp 


P(0) = 1, 


d( 


( 0 ) = 0 . 


A simple expression which satisfies equation ( |29| ) exactly, but violates the boundary conditions in equations (p 
used in models due to its simplicity, is that of a singular isothermal sphere, for which 

p(r) = 


and 


2t rGr 2 ’ 


1 GM{r) 


— constant. 


(27) 

(28) 

(29) 

(30) 
often 

(31) 

(32) 


However, any solution of both equations (|2£|) and (|30| ) must have a core (i.e . ro ^ 0), and, while equations (|ll|) and (|32[) 
correspond to the asymptotic limit of the solutions of equations (^9|) and as ro —> 0, they are not good approximations 
in general. We will compare our results for the nonsingular case with that of the singular isothermal sphere in what follows. 

The isothermal sphere has infinite mass (M oc r for large radii). Hence, in order to describe some realistic finite structure 
in terms of this model, we must truncate the sphere at some radius, rt . The total mass Ado of the isothermal sphere is then 
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M 0 = M(r t ) = / 4n p{r)r 2 dr = 4nporoM(g t ) 


where (t = rt/ro and M is the dimensionless mass: 


M(« = 


m <«■ 


4vrrgpo 

Using equation (|2q), we obtain 


M(r t ) = 


r 


G(47rGp 0 ) 1/2 

where Mt = M(C,t)- Solving for a 2 , we obtain 


while equation ( |26| ) yields the gas temperature: 
^ \ Vs 

t = ^W) 1/3 

kb 


MoV 

Mt/ 


(33) 

(34) 

(35) 

(36) 

(37) 


The existence of a truncation radius leads to the necessity of an external pressure to keep the system in equilibrium. This 
requires a significant correction to the form of the virial theorem stated in the previous section and, accordingly, to all of the 
other results obtained there. 


4 THE VIRIAL THEOREM FOR TRUNCATED ISOTHERMAL SPHERES AND ITS 
CONSEQUENCES 

Let us consider the general isothermal sphere density profile p(r), obtained in the previous section. From the ideal gas law, 
the pressure inside as a function of the radius is 


P(r) = ——p{ r ) = <? 2 p(r), 
m 

and at the outer edge 
pt = p(rt) = (J 2 p{r t ). 

The mean pressure inside the sphere is 

- fpdV 3/ 0 Ct p(CK 2 dC 3 M(Ct) _ 

P = r — = - n - Pt = ToTTTTPt = «(C t)Pt 


fdV 


c 4 3 p(c*: 


(?p{CtY 


(38) 

(39) 

(40) 


where = ry/ro. A jilot of a(Ci) is shown on Figure gj. The limit = oo, a = 3 is the case of a singular isothermal sphere, 
which will be considered in more detail below. 

The virial theorem for a static sphere in the presence of a surface pressure reads 


0 = 2K + W + S„ 


(41) 


where W is the gravitational potential energy, K is just Uth, the thermal energy, and S p is the surface pressure term. The 
thermal energy for a gas with ratio of specific heats 7 = 5/3 is given by 


Uth = \ J pdV= | a (Ct)ptV, 

where V is the total volume. The surface term is equal to 
S p = — / pr ■ dS = —3Vp t - 


Hence the potential energy is 

W = -2 U th ~S P = -2 ?—±Ut h 
a 

and the total energy is 
Em?—^U th . 


(42) 


(43) 


(44) 


(45) 
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Figure 1 . The ratio of the mean pressure to the boundary pressure a = p/pt for an isothermal sphere truncated at radius rt (solid 
curve), is plotted against the dimensionless truncation radius £ t = rt/ro. The ratio of the temperature of the isothermal sphere to that 
of a uniform sphere with the same mass and total energy according to the SUS approximation, b t = T/TguSi 1S plotted as well (dashed 
curve). 


In order to satisfy the virial theorem, the condition that Uth > 0 (and E < 0) requires a > 2, while Figure |l] indicates that 
a < 3.7846 in general. The virial temperature of the isothermal sphere is 

2 a GM 0 m 

T - S(c-2 ) fa,r„ ' (46) 

which we shall henceforth refer to as Tt/sQ- Since a/(a — 2) > 1 for any a, the temperature Ttis is always higher than 
Tsus, the standard value for a uniform sphere shown in equation (J^). If we write Ttis = brTsus, we find br > 2.13. This 
temperature correction factor is also plotted in Figure ^]. 

The virial radius r v n in this case is just the size of the truncated isothermal sphere, or r v i r = rt- For comparisons with 
the results of dynamical calculations of the formation of such an equilibrium object, we should interpret r v i r (and rt.) as the 
radius inside which hydrostatic equilibrium holds ( Cole fo Lacey 199t| ). Care should be taken not to confuse this definition of 
r v ir with the one that is most often used previously when identifying collapsed haloes in simulations, that of the radius inside 
which the mean overdensity is 18rr 2 « 178, as predicted by the top-hat solution if r v i r = r m /2. We will demonstrate below 
that our solution for r v i r = rt is not far from this value, even though the structure of the TIS is very different from that of 


the SUS approximation. This is consistent with results obtained using numerical simulations (Cole & Lacey 1996). 


For some purposes, the approximation of a truncated, singular isothermal sphere (“SIS”) has been used, instead of the 
exact solution, since it is algebraically simpler and the results can be obtained analytically, even though this approximation 
cannot be correct at the origin. For comparison with the nonsingular TIS results and as further illustration of the importance 
of the correction for finite boundary pressure discussed above, we calculate the virial temperature for such truncated, but 
singular isothermal spheres as follows. Using equation (Bl|) and the ideal gas law, the pressure at radius r is given by 


p(r) = p t 


( 7 ) 


where the boundary pressure at r = rt is given by 


pt 


k B T 


-pt 


2 

a pt = 


2-KGr? ' 


The average pressure inside the truncated isothermal sphere is then just 
p = 3 pt, 

(i.e. a = 3). The thermal energy in this case is given by 


(47) 


(48) 


(49) 


* Henceforth the notation “TIS” shall refer to a solution of the isothermal Lane-Emden equation with nonsingular boundary conditions 
of equation (BOI) at the origin 
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rr 3 [ 9 3 a\t 3 GMq 

Uth = -j pdV = - Pt V = — = -— 

while the surface pressure term in the virial theorem equation is given by 

2cr 4 rt _ 1 GMo 

~G _ ~2 rt 

According to the virial theorem, the potential energy is, therefore, given by 

W = - 2 U th -S p = -^r- = -\u th = 

G 3 r t 

and the total energy becomes 


E = U th + W = -^ = - 


1 GMq 
4 


( 50 ) 


(51) 


(52) 


(53) 


3 4 r t 

The SIS virial temperature is calculated as above for the nonsingular case according to equation (phi), except that a = 3 for 
the singular case, which yields 

_ 6 GM 0 m 

SIS 5 k B r m ’ ( ^ 

This virial temperature is a factor of bx = a/{a — 2) = 3 times larger than Tsus in equation © given by the SUS 
approximation. Equations ( pill ) and ( [5j| ) then yield the truncation radius for this SUS case, rt = (5/12)r m . Hence, in the 
limiting case of a singular, truncated isothermal sphere, the correction to the virial temperature resulting from the finite 
boundary pressure term is of considerable importance, while the difference between the actual radius of the isothermal sphere 
and the virial radius of the SUS approximation (i.e. r v ir,sus = r m /2) is more modest. In what follows, we shall derive the 
corrections for the actual, nonsingular truncated isothermal sphere which results from the top-hat collapse. 


5 CHOOSING A UNIQUE PROFILE: THE MINIMUM-ENERGY SOLUTION 
5.1 The solution 


As shown above, the truncated isothermal sphere solutions form a one-parameter family, described by (t = rt/ro - the 
truncation radius in units of the core radius. Specifying (t, the total mass, and total energy completely determines the solution. 
Alternatively, we can specify the mass, the total energy and the applied external pressure pt . Using these parameters, we can 
implicitly express the truncation radius as follows. The dimensionless density at is 
mpt 


p(Ct) 

- El - 

S 0 

II 

po 

V po 

From 

equations (|2g|^ 


Mo 

mG 

ro 

Mt 

k B T’ 


1 

( m* y 

Po 

47T 

\Mo 


k B T po 


(k B Ty 
V mG ) 


(55) 


(56) 


and, by substituting this po in equation (b5|), we obtain 




G 


\k B Tj 

Using equation (|45|), the total energy can be written as 
E= 2-a3k B T Mo 


a 


2m 


and, together with equation (|57|) , we finally obtain 

3 
2 

where 


E = ^ (4nG 3 p t M^ 1/4 E(( t ), 
2 - a(C t) 


m) = 


1 1/4 ‘ 


«(Ct) [p((t)M 2 ((t)\ 

Similarly, we can define a dimensionless external pressure p\ 


(57) 


(58) 


(59) 


(60) 
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Figure 2. Logarithms of the dimensionless energy E{C,t) (solid curve), external pressure p(C,t) (long-dashed curve), and radius |A^;(^t)| 
(long-dashed - short-dashed curve) versus dimensionless truncation radius for a truncated isothermal sphere. The vertical dot-dashed 
line indicates the value = 4.738 below which E > 0 and a TIS cannot exist. The vertical dotted line indicates the value = 2 9.4 at 
which the energy is minimal. The vertical thick line indicates the boundary of stability for the TIS solution (ft = 34.2, see § 5.2). TIS 
solutions are unstable to the right of this line. 


p(C*) = [£(C0] 

where 

(2£/3) 4 


Pt = 


4ttG 3 M 0 6 




(61) 


(62) 


We have plotted the dependence of the dimensionless total energy E and the dimensionless external pressure p on ft in Figure 
In order to indicate the dependence of the size of the sphere on ft, we nondimensionalize the radius rt, according to 
rtE _ 3 a — 2 ft 

2 _ 


— A e = — 


GMq 


M t 


(63) 


and plot —As against ft in Figure g, as well. This definition of As also corresponds to the familiar dimensionless energy 
parameter used in discussions of the stability of isothermal spheres (see § below). 

For a given mass Mo and external pressure pt, the solution is specified uniquely only if we can uniquely identify a special 
value of ft, or equivalently, of E. Apparently, for any truncated isothermal sphere of mass Mq which is confined by a given 
external pressure pt, there is a unique value of f t which minimizes the total energy E. We shall make the reasonable anzatz 
that this minimum-energy solution is the unique TIS solution preferred in nature as the outcome of the virialization of the 
sphere in the presence of a fixed external pressure. We offer evidence to support this anzatz in § [^for the cosmological top-hat 
problem. It is possible that this is a general result for any TIS formed by relaxation in the presence of a fixed external pressure, 
but we are only concerned here with spheres that evolve from cosmological initial conditions. 

The minimum value of E as a function of ft for a given pt is found by minimizing the dimensionless energy 75(ft) in 
equation (lso|). This occurs for ft = 29.4, for which E((t) = —0.2816, a(ft) = 3.73, Mt = 61.485, and p(ft) = 1.946 x 10~ 3 . 
For a given mass Mo and boundary pressure pt, the full details of the TIS solution are described in terms of these quantities 
as follows. The total energy E is calculated from equation (^), resulting in 

\(rTr){-^p) G 3/4 Pt 1/4 M 0 3/2 =-0.795 G 3/ V 4 M 0 3/2 . (64) 


Em 


The temperature T is obtained from equation ( J58| ) and (p£>[) , according to 

1/4 

/ 4-7T 1 

T 




m ^/ 4 


mT = 14.4 ^G 3/4 Pt 1/4 


M, 


1/2 


\p t M 1 J ks rt u 

Finally, the core radius ro and the central density po are calculated from equations (|5' 

1/4 


and (pa) 


r o 


\ 471 M t 2 


G 1/4 p - 1/4 M i/2 = 142 x 10 —2 g 1/4 p7 1/4 M ( 1/ 2 


0 


(65) 


( 66 ) 
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and 


P o 



1/4 

g- 3 /4 p? /4 m- 1/2 


449.5 G~ 3/4 p 3/4 M" 1/2 . 


( 67 ) 


The size of the TIS is then just 


r t = (tro = 29.4 r 0 (68) 

with ro given by equation (^6|). The dimensionless radius for this value of (jt corresponds then to Xe = —0.3326. 

Apparently, the core of the truncated isothermal sphere is very small compared to its overall size. Nevertheless, the profile 
is quite different from that of a singular isothermal sphere, in that it has a higher average pressure, lower temperature and a 
core. 


5.2 Stability 


The problem of the stability of truncated isothermal spheres has been discussed by several authors (e.g. Antonov 1962, 
Lynden-Bell and Wood 1968, Katz 1978, Binney and Tremaine 1987, Padmanabhan 1989, 1990). According to this literature, 
our unique minimum energy solution above is stable, as follows. The case at hand is that of a sphere bounded by a thermally 
insulating wall, distinct from the case of a sphere in contact with a heat bath. The family of allowed TIS solutions of the 
Lane-Emden equation can be parameterized by one variable, or, equivalently, by the combination of parameters given by the 
ratio of densities measured at ro and rt, TZ = po/pt and the dimensionless energy A e = rtE/GM 2 . For gravitationally bound 
spheres (E < 0), there are no equilibrium solutions for A e < —0.335. Only for A e > —0.335 are isothermal spheres possible. 
For large enough density contrast between the center and edge of the sphere, however, equilibria exist but are unstable. In 
particular, as first discovered by Antonov (1962) and later discussed in detail in Lynden-Bell and Wood (1968) - who called 


this “the gravothermal catastrophe” - isothermal spheres with A e > —0.335 and TZ > 708.61 (corresponding to (t = 34.2) can 
exist as solutions of the Lane-Emden equation but are unstable. Such spheres are metastable configurations, for which entropy 
extrema exist but are not local maxima. Only in the intermediate regime, where A e > —0.335 and 7 Z < 708.61, can isothermal 
spheres both exist and be stable. In this regime, the TIS can be a state of local maximum of the entropy. We note that our 
unique, minimum-energy solution for the TIS in the presence of a fixed boundary pressure pt, described above, satisfies the 
stability requirement for avoiding the gravothermal catastrophe, since (t = 29.40, TZ = [p(£t)] _1 = 514, and A e = —0.3326. 

Finally, we note that in the TIS solution derived here, the gas temperature is spatially uniform but the equation of 
state of this gas is not an isothermal one (i.e. changes of density for a given fluid element do not leave the temperature 
unchanged). In the latter case, the stability is analyzed by considering isothermal spheres in contact with a heat bath, rather 
than the thermally insulated spheres described above. For a gas with an isothermal equation of state, isothermal spheres are 
found to be unstable for density contrasts below that of the gravothermal catastrophe limit. This analysis does not apply to 
our minimum-energy solution above, however. Hence, we do not need to consider the classical results of Ebert (1955, 1957), 
Bonnor (1956) and McCrea (1957) which deal with the problem of gas spheres bounded by a given surface pressure in thermal 
equilibrium with an external heat bath, which differs from the present case. 


6 THE TRUNCATED ISOTHERMAL SPHERE FORMED BY TOP-HAT COLLAPSE 

The unique TIS solution found in § by minimizing the total energy of the sphere for a given boundary pressure pt can be 
applied to the top-hat problem, as follows. The virialized object which results from the collapse of a given top-hat density 
perturbation must have the mass of the top-hat and the same total energy as the top-hat before it collapsed and virialized. 
Fixing the mass Mo and the energy E of the TIS which results from the collapse of the top-hat amounts to fixing the boundary 
pressure pt as that value which makes the TIS with this total energy E correspond to the minimum-energy solution described 
above. Hence, the TIS which results from top-hat collapse must have = 29.4 and a = 3.73, as found in § This means 
that the virial temperature Ttis is larger than that found by the SUS approach, T$us, by the factor br = a/(a — 2) = 2.16. 
A plot of p/psus as a function of r/rt, demonstrating the significant difference between the two density profiles is shown on 
Figure || 

This conservation of the energy E of the top-hat before and after its collapse and virialization assumes that there is no 
extra pdU-work that needs to be taken into account due to the presence of the external boundary pressure pt which would 
otherwise alter the final total energy E compared with its initial value before collapse. This is appropriate for the case at hand, 
since the collapse prior to the epoch of virialization is that of a cold, pressure-free gas which collapses supersonically. As in the 
well-known, self-similar, spherical infall solution of Bcrtschinger (1985), the original energy is converted from potential energy 
at maximum expansion to a mixture of infall kinetic energy and potential energy during infall, with a negligible share of the 
energy going into compressional heating. In that solution, the infall is halted by a strong shock. At this shock, the kinetic energy 
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Figure 3. (upper panel) The postcollapse equilibrium density profile of a spherical top-hat density perturbation versus the radius in 
units of the radius of the tophat at maximum expansion for the minimum-energy TIS (solid line), the singular isothermal sphere (SIS, 
dashed line) and the uniform sphere approximation (SUS, dotted line). All densities are plotted in units of I Hit 2 p(t rn \\), the postcollapse 
density of the SUS approximation, (lower panel) The logarithmic slope of our minimum-energy TIS profile versus radius as compared to 
the constant slope of —2 for the SIS profile. 


of infall is converted primarily into the thermal energy of the shock-compressed gas, with only a small portion remaining as 
kinetic energy of subsonic, postshock infall. By analogy with this infall solution, therefore, we identify the boundary pressure 
pt in the case of our TIS solution, not as a fixed external pressure which acts on the top-hat boundary throughout its collapse 
and virialization, but rather as something like the instantaneous post-shock pressure in the infall solution. As such, it has a 
physical origin in the conversion of the original energy of the collapsing top-hat, itself, from potential energy at maximum 
expansion into kinetic energy of infall during collapse and, finally, into thermal energy of the post-shock gas, always conserving 
the original energy E of the top-hat. 

The final size of the sphere after collapse and virialization is related to the size of the parent top-hat density fluctuation 
at the epoch of turnaround by a collapse factor p defined by 


V=— ■ (69) 

T m 

Using equations (^3{j) and (^), it is easy to see that 


= _5 5 q(CQ - 2 Ct 

77 3 B 2 a(Ct) M(Ct) 


(70) 


A plot of ??(Ct) versus is shown in Figure ^j. For the minimum-energy solution, using f = 29.4 in equation JtoI ) yields 
p = 0.554, i.e. ry = 0.554 r m . This value of p derived for the minimum-energy TIS solution above is similar to the value usually 
cited for the postcollapse, virialized object in the SUS approximation, for which r v i r = r m / 2, corresponding to psus = 0.5. It 
is somewhat larger than the value for the limiting case of a singular isothermal sphere (“SIS”), for which p = 0.4167. 

The radius rt of the TIS solution is encountered by the outer boundary of the collapsing top-hat at a time somewhat 
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Figure 4. The collapse factor relating the radius of the initial top-hat at its point of maximum expansion to its radius after collapse 
and virialization, ■// = rt/r m , as function of the truncation radius The dotted line indicates the value of // given by the self-similar 
infall solution (see :j^i. 


earlier than the time t C oii of infinite density in the top-hat solution. We shall refer to the time at which the radius of the 
collapsing top-hat equals rt as t c ross, to distinguish it from t co u- At time fcross, the top-hat exact solution ([]]) yields an 
overdensity of 5 = 5 cross ~ 100, while the extrapolated linear solution at this time yields Sl = 5l, cross ~ 1.52. If we assume 
that the TIS forms instantaneously at t = fcross, then its mean density corresponds to a mean overdensity <5 = <5 cr oss ~ 100 
when compared with the background density at Across, and corresponds to a mean overdensity S = <5 C oii ~ 130 when compared 
with the background density at z C oii- This latter value differs somewhat from the conventional value of 187r 2 ~ 178, which is 
found in the SUS approximation. 

In applications of the top-hat model involving the Press-Schechter approximation, it is customary to identify the charac¬ 
teristic time of formation of objects of a given mass in terms of S c in equation (j|). Since <5 l, cross < S c , it may be appropriate 
to replace S c in such applications by the value 8l, cross ~ 1.52, implying that the number of objects formed at any epoch with 
a mass greater than some value will be somewhat higher for the TIS solution than previously assumed in applications of the 
Press-Schechter approximation. We note that comparisons of the predictions of the Press-Schechter approximation for the 
mass function of virialized haloes at any epoch with the results of IV-body simulations of cosmic structure formation from 
scale-free initial conditions have sometimes indicated better agreement if the standard value of 5 C = 1.69 is replaced by a 
lower value of 5 C (e.g. Efstathiou and Rees 1988; Carlberg and Couchman 1989; Klypin et al. 1995; Bond and Myers 1996; 
Cen 1998). Perhaps our TIS result that 5l, cross ~ 1-52 is part of the explanation for this. 

In order to check our derivation of the unique TIS solution described above and confirm the validity of the minimum- 
energy argument on which it is based, we take a second approach to the same problem, in the following section. In this second 
approach we match the properties of the final TIS equilibrium configuration to those of the post-shock gas in the well-known, 
frequently applied, self-similar, cosmological, spherical infall solution (Bertschinger 1981). As we shall see, this approach not 
only confirms our derivation above, but also provides a dynamically self-consistent origin for the final equilibrium state of the 
TIS. 


7 APPLICATION OF THE SELF-SIMILAR SPHERICAL INFALL SOLUTION 

The description of the postcollapse virial equilibrium of a top-hat density perturbation as a TIS with the same total mass 
and energy as the parent top-hat is not unique without the use of additional information which establishes the radius, the 
temperature, or the boundary pressure. We have argued above that a natural outcome of the dynamical relaxation of the 
collapsed top-hat is that unique TIS which makes its energy E correspond to a minimum-energy solution at fixed boundary 
pressure. This, in turn, establishes what the value of pt is for the TIS which results from the collapse of a given top-hat, and, 
thereby, fixes the unique TIS solution. To test this anzatz, we shall now place our top-hat density perturbation and its collapse 
and virialization in a dynamically self-consistent context, by comparing with the closely related problem of the spherical infall 
which results when a point-mass perturbation is added to a matter-dominated Einstein-de Sitter universe with initially 
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uniform density, for which an exact solution is possible. This latter problem represents the simplest kind of generalization 
of a uniform top-hat initial condition which can lead to an interruption of the collapse of the entire perturbation to infinite 
density at a finite time, due to the radial nonuniformity of the initial conditions. As such, this problem can be used to provide 
some insight into how inhomogeneities in a collapsing top-hat might halt its collapse prior to its reaching infinite density, 
causing it to thermalize the kinetic energy of collapse and relax to an equilibrium configuration. The spherical symmetry of 
this point-mass perturbation problem, moreover, means that one can identify a corresponding top-hat perturbation for the 
mass interior to any spherical shell in the point-mass perturbation problem, by defining the top-hat which contains the same 
total mass and energy inside the same size sphere at the same initial moment. Hence, a quantitative comparison of the solution 
of this problem and a suitably chosen, matching top-hat problem should be possible. We describe this in what follows. 


7.1 The self-similar spherical infall solution 


Consider the self-similar spherical infall solution for the flow which results when a point-mass perturbation is introduced at 


some point in an otherwise unperturbed Einstein-de Sitter universe (Bertschinger 1985). For a collisional fluid, this solution 
involves an accretion shock which occurs at radius rs = Xsrta, where rta is the radius of that spherical shell of matter 
surrounding the initial perturbation which is just turning around at any epoch t and where As is a constant. For an ideal 
gas with ratio of specific heats 7 = 5/3, the constant As ~ 0.34. As we shall describe in more detail below, it is appropriate 
for us approximately to identify the region of this self-similar infall which is bounded by the accretion shock at radius rs 
with our postcollapse equilibrium sphere. In order to do this, we must make the appropriate correspondence between the 
top-hat perturbation in our initial conditions and the “matching” top-hat in the infall solution which results when we add 
the initial point-mass perturbation to the unperturbed gas in a finite sphere so as to encompass the same mass in the infall 
solution as in our initial top-hat. The correct “matching” of the two solutions must ensure that the spherical shell which is 
just reaching the shock in the infall solution at some time t was following the same trajectory as the outer boundary of our 
top-hat perturbation at all times prior to the time it encounters the shock at rs- We outline this “matching” solution below. 
We argue, further, that this collisional fluid solution should apply equally well to the case of a TIS in a gas of cold, collisionless 
dark matter or any combination of fluid with 7 = 5/3 and cold, collisionless gas, as follows. Interestingly enough, for the 
same kind of perturbation in an unperturbed, purely collisionless, initially cold gas (e.g. cold dark matter), the self-similar 
solutions of Bertschinger (1985) and Fillmore and Goldreich (1984) yield an inner region of shell crossings, bounded by the 
outermost caustic located at r c , 1 = A c ,r rta, where A c ,i « 0.36. As a result, the outer boundaries of the shocked gas and of 
the region of the shell-crossings and caustics in the collisionless dark matter are almost coincident, and their respective mas s 
distributions are unbiased with respect to each other. In that case, we can focus in what follows on the Bertschinger ( 1981 ) 
solution for a collisional fluid with 7 = 5/3 and fi = 1, but bear in mind that our final result for the TIS applies equally well 
to a collisionless dark matter halo or to a combination of the two, as well, as long as fl = 1 . 

Following Bertschinger (1985), we define the dimensionless variable A = r/rt a (t), where rta is the radius of the shell 
that is just turning around at the time t. The fluid variables pressure p, mass density p, and mass m(r) interior to radius 
r, are non-dimensionalized using the following definitions of dimensionless pressure P, density D and mass M interior to 
dimensionless radius A = r / rta- 


p(r,t) = pb(j-j-^j P{ A), p{r,t) = p b D( A), 


m(r,t) = ^ Tvp b rt a M (\), 


(71) 


where p b is the mean, unperturbed background density at time t. Inserting these definitions into the fluid conservation partial 
differential equations leads to a set of ordinary differential equations. These equations, together with strong, adiabatic shock 
jump boundary conditions, the exact, pressure-free infall solution for the flow exterior to the shock and the assumption that 
the energy of the gas bounded by the shock is the same as the initial energy of the same mass when it was outside the shock, 
describe a self-similar infall solution with an accretion shock. Bertschinger (1985) presented detailed results of his numerical 
solution of these equations both graphically and in tabular form. However, since we shall need a finer grid of solution values 
than is tabulated there, we have re-integrated the resulting equations ourselves for the current paper, using the stiff equation 


solver RADAU5 by E. Hairer and G. Wanner (Hairer, Norsett k, Wanner 1996). 


The shock occurs at fixed A = A s = rs/r ta and propagates according to 
rs{t) = Xsr ta (t) oc t s/9 . 

The mass inside the shock is determined by: 
m(r s ,t) = ^tt p h rta{t)M (\ s ), 

The self-similar infall solution yields, for 7 = 5/3: 

A s = 0.338976, M (As) = 3.78759, P(As) = 9.73563, 


D(Xs) = 16.5341. 


(72) 


(73) 


(74) 


© 0000 RAS, MNRAS 000, 000-000 






14 


P.R. Shapiro, I.T. Iliev and A.C. Raga 



Figure 5. The ratio of pressure force to gravitational force inside the accretion shock in the self-similar infall solution. A value f = 1 
implies hydrostatic equilibrium. The shock is located at log 10 A $ = —0.469831. 


The mean overdensity inside the shock relative to the mean unperturbed background density is constant in time, and has a 
value 

M(\ s ) 


m(r s ,t) 

PS 4 o 

3 nr sPb 


A| 


= 97.24. 


(75) 


Consider the post-shock gas in this self-similar infall solution. Strictly speaking, this postshock flow is neither isothermal 
nor in hydrostatic equilibrium. As we shall see, however, the variation of temperature with radius is rather weak, so isother- 
mality is a reasonable approximation and the post-shock gas is in approximate hydrostatic equilibrium, as well. We shall 
demonstrate this as follows. If the gas is in hydrostatic equilibrium, it must satisfy the equation 


1 dp 

—— = -V< i>. 
p dr 

Let us define a dimensionless ratio 


(76) 


(77) 


which measures the degree to which hydrostatic equilibrium is obtained at any radius r, so that f = 1 if equilibrium occurs. 
In order to demonstrate how close the exact solution is to hydrostatic equilibrium inside the shock, we rewrite § in terms of 
the dimensionless fluid variables, according to 


C(A) 


9 [ A 2 P'(A) 

2 M(X)D(X) 


(78) 


We use the exact numerical solution to plot £(A) in Figure Apparently, the postshock gas is always close to hydrostatic, 
with a maximum departure of £ from unity of 9 per cent just behind the shock, where the gas still has a small infall velocity 
and the pressure is a little lower than the value required to balance gravity. 

The approximation of equilibrium becomes perfect as we approach the origin. We can demonstrate this analytically, using 
the simple asymptotic scalings of the fluid variables with distance there. In terms of A, the solution obeys the following scalings 
at small A: 


P = p( A)A _5/2 , D = D{ A)A“ 9/4 , M = M( A)A 3/4 , (79) 

where P, D and M are only weakly dependent on A for small A, P(0 ) = 2.40676, D(0) = 2.60173, and M(0 ) = 10.4069. 
Inserting equation (|t^) into equation ((F^) for £(A) shows that £ is independent of A as A approaches 0, with £(0) = 1. 

The postshock gas is also not far from isothermal, despite the strong increase of pressure and density towards the centre. 
The volume-averaged temperature interior to the shock radius is related to the immediate postshock temperature Ts according 
to 


(T)v 

T S 



P(X)/P(\s) 

D(X)/D(X S ) 


X 2 dX. 


(80) 
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According to the scaling laws in equation (|?^), valid at small radii and approximately correct throughout much of the postshock 
region, temperature varies with radius only as r -1 ^ 4 . In terms of the scaling laws, 


(T)v „ Ha -!/ 4 p(Q)/p(*s) 

Ts ~ 11 s D(0)/D(\ s ) 


(81) 


while the exact solution yields 


(T)v 

T s 


1.386, 


(82) 


which is remarkably close to unity. 

In short, the self-similar infall solution indicates that the dynamically self-consistent evolution of an uncompensated, 
spherically-symmetric, positive, cosmological density perturbation in a cold gas leads to continuous infall and the creation of 
a continuous sequence of spheres of collapsed gas which are, at any instant, in a state which is close to that of an isothermal 
sphere in hydrostatic equilibrium, bounded by the accretion shock. In what follows, we identify the particular self-similar, 
shock-bounded sphere in this solution which instantaneously corresponds to the postcollapse, virialized sphere which forms 
as a result of the equivalent top-hat perturbation. 


7.2 Matching the post-shock gas in the infall solution with the postcollapse TIS of the top-hat 


In what follows, we attempt to find the TIS solution for the postcollapse virial equilibrium of a top-hat density perturbation 
which most closely matches the characteristics of the accretion-shock-bounded sphere in the self-similar spherical infall solution 
when the mass and energy inside the latter sphere equal the mass and energy of the top-hat. We require that the trajectory 
of the top-hat outer boundary be the same as that of the immediate postshock fluid element in the infall solution during the 
preshock infall of the latter. This fixes the radius of the shock and the time at which the shock encloses a mass equal to that 
of the corresponding top-hat. By comparing this shock radius with the radius of the immediate postshock fluid element at its 
point of maximum expansion, we shall derive the effective collapse factor for the corresponding top-hat from its maximum 
expansion to its postcollapse virial equilibrium state, by identifying the shock radius rs as the radius rt of the matching TIS. 
As we shall see, this infall shock radius in almost precisely equal to the radius derived previously for the minimum-energy TIS 
solution, although it is necessary to make a tiny downward adjustment of the value of rt relative to rs in order to match a 
TIS solution with the same total mass and energy. In particular, an isothermal sphere with exactly the same mass and energy 
as this shock-bounded sphere in the infall solution can only be in hydrostatic equilibrium (i.e. a solution of the Lane-Emden 
equation) if its radius is a tiny bit smaller than the shock radius in the infall solution. We shall make an approximation, 
therefore, to determine the actual radius rt, slightly below rs, of the equilibrium TIS solution which most closely matches 
the shock-bounded sphere in the infall solution, as follows. As we shall see, the resulting TIS is not far from our TIS solution 
derived above from the minimum-energy argument, so we can be confident that our minimum-energy TIS solution has a sound 
physical basis. We begin by solving for the shock radius rs and time t cr oss at which the collapsing outer boundary of the 
top-hat reaches the radius of the shock in the infall solution, when the gas enclosed by the shock in that solution contains the 
same mass Mo and energy E as the top-hat. 

The total mass of the top-hat we consider is: 

Apr q 

M 0 = — r m(l + S m )p b (tm) ■ (83) 

We require that the masses of the top-hat and the post-shock gas in the self-similar infall solution be equal at the time when 
the outer radius of the top-hat crosses the radius of the shock. This yields 
4 3 

Mo — rri (rs ■ Across) — “ 71 pb (fcross) r t a (fcross) Af (A<? ) , (84) 

and, therefore, 


rta (fcross) — 


M 0 


1 1/3 


M(As) ^p b {t C 


(85) 


where 


pbiT) — (67lGf ) , Pi> (fcross) — Ph()( 1 "T Across) , 

and pbo = 3 Hq /8nG, the mean background density at present. Combining equations (| 


, we can write 


r ta {t cross) 


3 (36^)!/ 3 G 2/3 

2 [M(As)] 2 / 3 (i+Zc 


s )p, 


1/3 


and 


( 86 ) 


(87) 
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( 3 A s jfl/3/-, , -,-1 -1/3 

rm _ \47T/ ^[M(A S )]1/3 M 0 ( 1 + Zcross ) p »o ■ 

The radius of the shock wave at tcross is, using equations (|2|) and (|85]), given by 
r S (tcross ) = Xsrta (tcross) = C 1 + S rn ) 1/3 r m 

= r7 S M 1 n (! + <M V3 SL : cros3 r m = 0.362445i, cross r m , 
where 5l, cross is the linear density contrast corresponding to the time t c r0 ss- Expressed in terms of r m , this yields, 


( 88 ) 


(89) 


rs(t c 


Vsr„ 


(90) 


where r/s = 0.564393 (Bcrtschinger 1985). 

Our goal here was to use the shock radius rs(t c ross) to establish rt, the size of our postcollapse equilibrium TIS, by equating 
77 in equation ( |&s| ) with rjs in equation (|)(]). This value of rjs is extremely close to the value of 77 = 0.554 derived above for the 
unique minimum-energy TIS solution for the postcollapse equilibrium of the top-hat. Such an agreement is remarkable and 
can hardly be a coincidence. We take this as a strong confirmation of the validity of that minimum-energy solution. However, 
as mentioned earlier, we also note that this tiny difference between 77 = 0.564 and rj = 0.554 is nevertheless enough to put the 
isothermal sphere with mass M, energy E and radius rt = rs ever so slightly beyond the boundary of allowed solutions of the 
Lane-Emden equations, according to Figure^] and equation (|Fo|). In other words, since A e = — ( 3 / 5 )t 7 s = —0.338636 < —0.335 
in this case, such an isothermal sphere is just slightly too large to be an equilibrium TIS. In what follows we shall refine this 
estimate of the TIS radius r t so that it departs from r$ only minimally but nevertheless preserves the physical basis in terms 
of the self-similar infall solution. 


7.2.1 The TIS solution whose size is closest to the shock radius in the infall solution 

The significance of the remarkably close agreement between the radius rs (tcrosa) derived above and the radius rt of the 
minimum-energy TIS for the same parent top-hat is made clear by the curve of 77 versus in Figure ^ which shows that 
77 = 0.554, the value for (jt = 29.4, the minimum-energy TIS solution value is close to the maximum value of 77 for all values of 
£ t , while rjs = 0.564 is only just above the curve. Since the postshock flow is only approximately isothermal and hydrostatic 
it is reasonable to suppose that the final relaxation to equilibrium keeps the total energy and mass constant but allows the 
radius to shrink slightly, just enough to relax to a stable equilibrium TIS. The simplest way to refine the value of r t (or, 
equivalently, of 77 ) downward from rs (or rjs) so as to make an equilibrium TIS solution possible, therefore, is just to take 
77 = r/max , the maximum value of 77 versus t and, hence, the allowed value which is closest to r/s- This yields the estimate 
771 = r/max = 0.558, for which = £ t ,i = 34.2 and a = ai = 3.66. This rim ax solution corresponds exactly to the unique stable 
equilibrium TIS for which A e = —0.335 and TZ = 708.61. This value of 77 = 771 differs from rjs by only 1.3 per cent. With this 
choice of 771 , the TIS solution is extremely close to that from the minimum-energy solution, with 771 differing by only 0.7 per 
cent, a and T for the two solutions differing by less than 2 per cent for each, and (i.e. the core radius ro) by only 17 per 
cent. 

This agreement between our minimum-energy TIS solution and the infall accretion-shock size matching TIS solution 
confirms the validity of our minimum-energy TIS solution. We note, however, that even the very small difference found 
between the two TIS solutions is just enough to place the size-matching solution on the boundary, described above, of 
gravothermal catastrophe. Apparently, self-similar spherical infall leads naturally to the formation of a shock-bounded sphere 
which is close to that of an isothermal sphere on the edge of the gravothermal catastrophe regime. 


8 DISCUSSION AND SUMMARY 

8.1 Summary of the minimum-energy TIS solution 

We have argued above that a cosmological top-hat density perturbation in an Einstein-de Sitter universe which virializes 
following its collapse can be described self-consistently by a unique TIS solution which conserves the energy of the initial 
top-hat and is bounded by a pressure which makes that energy a minimum. The temperature derived in this way for the 
sphere in virial and hydrostatic equilibrium is larger by a factor of approximately two than the value previously derived 
by satisfying energy conservation and the virial theorem for a postcollapse sphere of uniform density for which the surface 
pressure term is neglected. The origin of the boundary pressure term in the TIS solutio n can be understood self-consistently 
by comparison with the well-known self-similar, spherical infall solution of Bertschinger( [l985l) for the problem of the collapse 
of matter onto a point-mass perturbation added to a uniform, cosmologically expanding background. In the latter problem, 
spherical mass shells surrounding the point mass decelerate relative to the mean cosmic expansion, reach a maximum radius 
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and collapse back toward the center supersonically, unaffected by pressure forces, following a trajectory identical to that of 
the outer boundary of a top-hat perturbation which encompasses the same total mass and energy. Unlike the collapse of the 
top-hat, however, the self-similar infall is interrupted by a strong accretion shock which decelerates the flow and thermalizes 
the kinetic energy of infall, allowing the postshock gas to settle subsonically into a state which is approximately hydrostatic 
and isothermal. By comparing the shock-bounded sphere in this self-similar infall solution to the unique minimum-energy 
TIS solution which we derived here by a completely independent line of reasoning, we have been able to confirm the validity 
of the TIS solution while explaining its dynamical origin and the nature of the boundary pressure. The boundary pressure, 
in short, is closely related to the presence of the accretion shock in the self-similar infall solution. While the TIS boundary 
pressure is not identical to the instantaneous postshock pressure in the infall solution (it is slightly higher), the TIS boundary 
pressure corresponds to the slightly adjusted postshock pressure required to transform the quasi-hydrostatic postshock flow 
into an exact hydrostatic equilibrium while conserving the total energy inside the sphere. We have also used the properties of 
the self-similar infall solution in the case in which both collisionless dark matter and a collisional fluid are present to argue 
that we can apply our TIS solution to either component or both. For the collisionless component, velocity dispersion takes 
the place of temperature. If both components are present, then the density of each component i, pi, is simply proportional to 
the total, pi = ptot^i, where fij is the universal density parameter for component i. 

The TIS solution for the postcollapse equilibrium of a top-hat of mass Mq which would collapse at epoch z co n if its collapse 
were not interrupted by virialization depends upon Mq, z co ii, h, the mean total background density pb,o (where pb,o oc h 2 ), 
and the ratio of truncation radius i't to core radius ?'o, according to 
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a — Z 

(91) 

Vc 

oc 

?? - i / 2 M i/ 3 (l + 2coii ) i / 2 /i i/3 ! 

(92) 

T m 

oc 

M^ 3 (l + 2 col i) - 1 h -2/3 , 

(93) 

Po 

oc 

— 3 "YY PbO ( 1 Zcoll ) ? 

T M t 

(94) 


where a and Mt depend only on (jt and where h is the Hubble constant in units of 100 km s -1 Mpc -1 . The radii rt = rjr m 
and ro = rt/Ct have the same scaling with Mo, 1 + z c oil, and h as does r m . 

The value we have derived above for £t for our minimum-energy TIS solution is C,t = 29.4. All of the parameters of this 
unique TIS solution are summarized in Table [j]. Using this value and evaluating the proportionality constants in the scaling 
laws (^J) (fTil), the scaling laws become the following: 
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The mean density inside the TIS is 

P = 3 ^ 7 T^Po = 7.25 x 10" 3 po = 18 tt 2 (^-p b o(l + Zcdi ) 3 = 130.5p w (l + z cM ) 3 ■ (102) 

St V 5 / St 

Together with a numerical solution of the Lane-Emden equation in dimensionless form for the sphere with parameters given 
in Tabic [l], shown in Figures and [], these scaling laws provide the complete description necessary to apply our TIS solution 
to some cosmological problem. 

For some purposes, an approximate fitting formula for the numerical solution of the Lane-Emden equation for isothermal 
spheres, equations (^) and ([l^), is convenient. Although we have worked exclusively with the actual numerical solution of 
equation (^) throughout this paper, so as to minimize our errors, we shall now provide such a fitting formula which is useful 
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Table 1 . Summary of the minimum-energy TIS solution 


Quantity 

Value 

Ct - a. 

29.40 

r 0 

II 

3.73 

M(Ct) . 

61.48 

b T = t tis 

T T SUS - 

2.16 

n - . 
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' rm 

r TIS 

1.11 

r SUS 
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1.52 


(Across). 
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^(^coll) . 

130.5 

Pt 

1.95 x 10“ 3 

PO 

77. 

514 

PO 

1.80 x 10 4 

Pfc>,coll 

Pt 

35.0 

Pb,c oil 


as a good approximation over a wider range of radii than has been presented in previous treatments in the literature. We 
adopt the form 


P _ A B 

Vo ~ a? + C 2 ~~ b 2 + C 2 ’ 


(103) 


with the parameters A, B, a , and b sensibly chosen (Natarajan and Lynden-Bell 1997; henceforth, NL97). In the Appendix, 
we briefly describe the method used to obtain such a solution, adjusted to the problem at hand. For the purposes of this 
paper a very good approximation is that of equation (103) with the following parameters: 


(A, a 2 ,B, b 2 ) TIS = (21.38, 9.08,19.81,14.62), 


(104) 


which differs from that in NL97 and is a better 6 t over a larger range of radii than theirs (see Appendix Al). It fits the exact 
solution within a fractional error of 3 per cent for 0 < £ < 40, which is a range that fully encompasses that of our final TIS 
solution, with a perfect match to the exact solution at the boundary g = g t = 29.4. 


8.2 Comparison of the minimum-energy TIS solution and the uniform sphere and singular isothermal 
sphere approximations 

The temperature derived here by the minimum-energy TIS solution for the postcollapse sphere in virial and hydrostatic 
equilibrium which follows top-hat collapse is a factor of approximately two larger than the value previously derived by 
satisfying energy conservation and the virial theorem for a postcollapse sphere of uniform density for which the surface 
pressure term is neglected, as summarized in § In particular, if we write Ttis = brTsus, then br = a/(a — 2) = 2.16. 
We shall address the question of the significance of this revision of the postcollapse virial temperature below. The size of the 
TIS sphere, given in terms of the radius r m of the top-hat at maximum expansion according to rt = 0.554r m , is actually 
not far from the size of the SUS sphere, r v i r = r m /2. However, in terms of the TIS core radius r o, the truncation radius is 
large, r t /ro = 29.4, implying that the TIS is very far from uniform density. Despite this relatively small core radius, the TIS 
solution is also quite different from that of a singular isothermal sphere; for the latter, the ratio of the average density to 
that at the surface is 3, while we find a value of 3.73 for the nonsingular TIS, and the truncation radius in the singular limit 
is only rt = (5/12)r m . Similarly, the correction factor for the virial temperature of the singular isothermal sphere relative 
to that of the uniform sphere is br = 3, as opposed to the value for the TIS which is br = 2.16. The central density of the 
TIS is 1.8 x 10 4 times larger than the mean density of the background at the collapse epoch for the initial top-hat. Finally, 
this solution predicts that the top-hat will virialize somewhat earlier than the nominal collapse time of the top-hat, since 
the outermost mass elements encounter a shock in the infall solution at finite radius. In terms of the usual overdensity of 
the top-hat predicted by linear theory, this implies a revised value of S c = ( 5 c )tis where ( 5 c )tis ~ 1.52, rather than the 
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standard value of 1.686 based upon extrapolating the linear growth to the epoch at which the nonlinear top-hat solution 
predicts infinite density. 


8.3 Comparison of the minimum-energy TIS solution density profile with results of numerical simulation 
of cosmic structure formation and with observations 

The hydrostatic isothermal sphere has been widely adopted as an equilibrium model for the density profiles of individual 
cosmic structures. Mass decomposition studies of galaxies based upon rotation curve data, for example, often adopt the 
density profile p oc [1 + (r/a) 2 ] -1 , an approximation to a nonsingular isothermal sphere, for the galactic halo (cf. Freeman 
1992). As another example, galaxy cluster X-ray brightness profiles are often interpreted in terms of the so-called “/1-model” 
in which the X-ray emitting gas is an isothermal sphere in hydrostatic equilibrium in the potential well of a virialized dark 
matter halo, for the purpose of measuring the cluster total mass density and baryon fraction profiles (e.g. Cavaliere and 
Fusco-Femiano 1976; David, Jones and Forman 1995). In this case, the gas density varies as p oc [1 + (r/a) 2 ]^ 3/3 ^ 2 , and the 
value of (3 = 2/3, which corresponds approximately to that of a nonsingular isothermal sphere, is consistent with the data on 
high temperature X-ray clusters (e.g. Jones & Foreman 1992; Arnaud and Evrard 1998). These examples of empirical fits are 
only loosely motivated by theoretical expectations for the dynamical evolution of the objects in question. 

Recently, Navarro, Frenk & White (1996, 1997; NFW) suggested, based upon the results of cosmological N-body simula¬ 
tions of the Cold Dark Matter (CDM) model, that dark matter haloes which condense out of the background can be fit by a 
universal density profile, of the following form: 

pNFW = 7 7-\"7T“- 7 -7X, (105) 

(r/rs)(l + r/rs) 2 

where the fitting parameters rs and ps are related to each other and reflect the different collapse epoch for different mass 
objects. This NFW profile shares the large-radius asymptotic shape, r -2 , of an isothermal sphere only for intermediate radii 
( r/rs > 1)) while steepening to r ~ 3 at large radii ( r/rs >> !)• At small radii ( r/rs « 1), the NFW profile flattens to r -1 . 
Unlike the nonsingular isothermal sphere, the NFW profile lacks a finite core. It is natural to ask, therefore, if the apparent 
differences between the TIS solution presented here for the postcollapse virial equilibrium of a top-hat density fluctuation and 
the NFW density profile are significant and, if so, what they mean. A full treatment of this question is outside of the scope 
of this paper. However, we will briefly address the issue as follows. 

There are several questions which must be answered in order to compare the TIS and NFW profiles in a meaningful way. 
First, are they solutions of the same initial value problem? Second, is the NFW halo a state of relaxed, virial equilibrium as the 
TIS solution is posited to be? Third, is the NFW profile really an accurate fit to the exact results for the dark matter haloes 
which result from the growth of density fluctuations in a CDM model universe. Fourth, is there a difference between the gas 
and dark matter profiles expected when objects condense out in the CDM model in general, or when the dark matter profile 
follows a NFW shape, in particular? Fifth, what do the observations of actual density profiles for cosmologically condensed, 
virialized structures say regarding the comparison between the TIS solution presented here and the NFW profile? 

Regarding the first question, the NFW profile is suggested to be the generic outcome of the condensation of haloes of 
collisionless dark matter in hierarchical clustering models like CDM, starting from Gaussian random noise density fluctuations. 
The TIS profile presented here is the outcome of the collapse of an uncompensated spherical top-hat, instead, albeit one in 
which the presence of some smaller-scale nonuniformity leads to relaxation to a final equilibrium state of finite density. We 
have argued that a solution very similar to our minimum energy TIS also arises if we consider the self-similar spherical infall 
solution of Bertschinger (1985), which is already close to isothermal and hydrostatic equilibrium inside the shocked region, 
and assume that rapid equilibration takes place within this shocked region. In that case, there is a continuous infall which 
leads to a continuous sequence of TIS solutions of ever larger mass as more mass encounters the accretion shock over time, 
and the shock represents a boundary which separates the relaxed region in equilibrium inside the shock from the unrelaxed, 
infalling matter at larger radii. As such, the TIS and NFW profiles do not necessarily arise from identical initial conditions. 

The simple spherical top-hat model has been quite successful, however, in characterizing the rate of formation of haloes 
of different mass at different times in hierarchical clustering models, as in the Press-Schechter approximation. The standard 
approach to this problem is to filter the density field in such models spatially with a spherical window function and to 
apply the top-hat model to predict the nonlinear collapse of initially linear-amplitude density fluctuations whose wavelength 
corresponds to the size of the window function. In a similar way, we might expect the TIS solution found here to be useful 
in characterizing the internal structure of the virialized objects which result from this nonlinear collapse. We note that there 
is currently no fundamental derivation or detailed theoretical explanation of the empirical NFW profile, so it is difficult to 
say just which initial conditions are necessary to produce it and how different these can be from that of the top-hat or of the 
spherical infall solution. 

Some attempt to address this was made recently by Huss, Jain and Steinmetz (1998), who added different amounts of 
velocity dispersion to the initial conditions of a spherically symmetric density perturbation evolved by N-body simulation to 
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see if the NFW profile would still fit the end result of its collapse. They found that the NFW was a decent fit for a range 
of initial velocity dispersions, but that it was possible to fit the simulation results better with a profile that was flatter at 
small radii than the NFW profile, as flat as r -0 ’ 5 . Syer and White (1998), on the other hand, have argued that the actual 
halo profiles in hierarchical clustering models are the product of the merging history of subclumps to form the halo, as an 
explanation for the cuspy inner density profile. They suggest that the slope of this inner density profile is different for haloes 
arising from linear density fluctuations whose initial power spectra have different slopes. If so, then the NFW is not a unique 
fit, since the actual profiles must depend upon the slope of the power spectrum at the wavelengths responsible for producing 
haloes of a given mass. This suggests that the haloes observed in N-body simulations are not fully relaxed, in the sense that 
they have not lost memory of the initial inhomogeneities which led to their formation. According to Syer and White (1998), 
in fact, the faster the initial density fluctuation power spectrum P(k) oc k n decreases with increasing k, the shallower is the 
predicted central density profile of the haloes which form on this scale, with p oc r r , and 7 « 3(3 + n)/(5 + n ). For the 
particular case of galactic scales in the standard CDM model, the effective power-law index is n « —2.5 and the scaling law 
would predict 7 sw 0.6, smaller than the NFW value. Only for larger-mass haloes, those at galaxy cluster scales, where n « —2, 
does this argument predict that 7 « 1 as in the NFW profile. We note that, according to the scaling law of Syer and White 
(1998), as n approaches —3, the approximate asymptotic value in the small wavelength limit in the CDM model, the inner 
profile is predicted to flatten to a constant core value. It is tempting to speculate that the scaling argument, if correct, might 
be pointing to our TIS solution for the virialized inner portion of a halo which condenses out of a CDM model universe, for 
smaller mass haloes, those which originate from density fluctuations in the small wavelength limit. In this limit, there are mass 
fluctuations of comparable amplitude on all scales, collapsing out nearly simultaneously, which may be the limit in which the 
phenomenon of repeated mergers of isolated, dense subclumps to build up a larger halo is not prevalent enough to generate 
a central density cusp. 

We have focused above on the differences between the NFW and TIS profile behaviour at small radii, but there is also a 
difference at large radii; the logarithmic slope of the TIS density profile drops as low as —2.5 at intermediate radii (r/r 0 ~ 9) 
and rises to —2.1 at the outer radius (see Fig®, whereas the outer slope of the NFW profile approaches —3 (i.e. at the outer 
radius, (r/rs)max, the NFW slope is 3 — 2/((r/rs) maz + 1)). This difference may reflect the transition between the inner 
virialized part of the halo, which we have characterized as isothermal and hydrostatic in formulating our TIS solution, and the 
outer region, involving continuous infall, which is neither isothermal, nor hydrostatic. The infall, itself, is generally predicted 
by spherical models to have a density profile closer to r~ 2 than to r -3 (e.g. Gunn and Gott 1972; Fillmore and Goldreicli 1984; 
Hoffman and Shaham 1985; Bertschinger 1985; Moutarde, et al. 1995). When account is taken of the presence of an accretion 
shock in the fluid case or of density caustics and shell-crossing in the case of collisionless matter, such as the self-similar 
infall solution of Bertschinger (1985), the density profile inside the shocked region and the region of shell-crossing is also 
closer to r~ 2 than to r -3 , even before any dynamical relaxation of the matter to an equilibrium profile occurs. Across the 
boundary between these two regions, the outer region of infalling matter and the inner region of shocked fluid and collisionless 
shell-crossing, however, there is a density jump, with a higher density inside the outermost density caustic than just outside 
this radius. To the extent that such a spherical model can represent the average behaviour of a more inhomogeneous initial 
condition, therefore, we expect the average density in the more general case to drop more steeply than r~ 2 across some 
transition zone at r ~ r v ; r where the actual transition in the collisionless matter is not discontinuous but is rather of finite 
radial extend. This would account qualitatively for the steepening in the NFW profile from r~ 2 to r^ 3 at large radii. 

The question of whether the haloes in the N-body simulations of the CDM model, fit by the NFW profile, are relaxed 
and in equilibrium has been partially addressed by Tormen, Bouchet and White (1998) for the case of collisionless matter in 
an Einstein-de Sitter universe in which the density fluctuations are Gaussian random noise with a power spectrum which is 
a scale-free power-law of index n = — 1. They report that the haloes, while not isothermal, are not very far from isothermal. 
They also find that the particle velocities at small radii are nearly isotropic while those at large radii are predominantly radial. 
As another approach, these authors compared different mass estimators for a set of haloes, all based upon an application of the 
Jeans equation for the dynamical equilibrium of a spherically symmetric collisionless system, allowing for different degrees of 
departure from the condition of velocity dispersion isothermality and isotropicity. They found that these mass estimators all 
compared favorably with the true mass if applied to estimate the total mass interior to the virial radius, while the estimator 
which took account of anisotropic velocity dispersion was in good agreement with the true mass at all radii. This suggests 
that the haloes are approximately and on average in a state of dynamical equilibrium but not fully relaxed to the point of 
having isotropic and isothermal velocity dispersion. As an additional probe, Tormen et al. (1997) also considered the fate of 
their simulated haloes at late times under the artificial circumstance in which a halo identified at some epoch is allowed to 
evolve further but with all matter removed outside of a sphere of radius twice the size of the original halo. This was intended 
to break the clustering hierarchy so as to allow the halo to evolve to a relaxed equilibrium state in the absence of further 
influx of perturbing substructure. Interestingly, the density profiles of the resulting haloes tended to flatten thereafter at very 
small radii, although the authors suggested that this might reflect some numerical error or artifact. Otherwise, the haloes 
showed less substructure with time but continued to be well fit, outside the central region, by the NFW profile. 

What do numerical N-body simulations tell us about the accuracy or universality of the NFW profile for the haloes that 
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form in the particular Gaussian random noise, hierarchical clustering model, CDM, which it was originally proposed to fit? 
Unfortunately, the results to date are not yet conclusive. For the purposes of discussion, it is useful to consider a broader 
family of density profiles, of the following form: 
po 
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(e.g. Hernquist 1990; Zhao 1996). This density profile approaches a power law of slope —7 at small radii, and a power law of 
slope — /3 at large radii, with a break at the characteristic radius rs and a transition region between whose width increases 
with increasing a. The NFW profile corresponds to the parameters (a, j3, 7 ) = (1, 3,1). An earlier suggestion of a galactic halo 
density profile for elliptical galaxies by Hernquist (1990) corresponds to (a,/3, 7 ) = (1,4,1). A co mpreh ensive review of the 
attempts to fit numerical simulation results for halo density profiles to those of the form in equation jlQC| ) above is beyond the 
scope of the present paper. We will mention only a few of the most recent and refer the reader to these papers and references 
therein for further description. Tormen et al. (1997) compare their simulation results by an N-body tree-code for haloes in 
the case of a scale-free power-law power spectrum of initial density fluctuations with n = — 1 to the NFW and Hernquist 
profiles, resolving individual haloes with 20, 000 simulation particles. Fitting the region between 0.01 and one times the virial 
radius (defined roughly in terms of the conventional SUS approximation for the average overdensity within a postcollapse, 
virialized top-hat), they find a somewhat better fit for the NFW profile than for the Hernquist profile. Recent high-resolution 
simulations of galaxy halo formation in the CDM model by Kravtsov, et al. (1998), however, by a new Adaptive Refinement 
Tree (ART) N-body code with a comparable number of particles per halo and a larger sample of haloes to analyze find a 
systematic deviation from the NFW profile at small radii, with a much flatter central profile, corresponding to p(r) oc r 1 with 
7 « 0.2 (i.e. (a, /3, 7 ) « (0.5, 3, 0.2); we note that Kravtsov, et al. define their own a as the inverse of the one we use here, so 
we transform their value to our notation). At small radii, at least, this is perhaps more consistent with the flat core of the TIS 
solution than with the cuspy center of the NFW density profile. On the other hand, very high resolution N-body simulations 
by a parallel N-body tree-code, reported by Moore, et al. (1998), of two rich-cluster-sized haloes in a CDM model in which a 
range of resolutions up to more than 10 6 particles within the virial radius was used, finds that the simulated profiles deviate 
at small radii from NFW in the opposite sense, becoming as steep as p(r) oc r 1 with 7 « 1.4 (i.e. (q,/ 3 , 7 ) « (0.7, 2.8, 1.4)), 
which they believe is reliably determined down to radii less than 0.01 times the virial radius. A qualitatively similar result 
regarding the density cusp is reported by Fukushige and Makino (1998) based upon their N-body simulations of the formation 
of a single, massive, galactic halo in the CDM model by direct summation of two-body gravitational forces utilizing the 
Grape-4 special-purpose computer. The full simulation volume of comoving radius 2 Mpc in present units contained 786,000 
particles and more than 10 12 Mq, but they do not indicate what fraction of these particles were contained within the virial 
radius of the final object. 

There are at least two possible explanations for the different results of these latest simulations. It may be that, as Syer 
and White (1998) suggest, the NFW profile is not universal but depends instead on the shape of the density fluctuation power 
spectrum at wavelengths responsible for the particular haloes in question. As such, the galactic halo profiles of Kravtsov, et 
al. are expected to be flatter at small radii than NFW, as described above, in contrast to the large, cluster-mass haloes of 
Moore, et al. (1998). This scaling argument alone does not, however, explain why the Moore et al. (1998) haloes are more 
cuspy even than those of Tormen et al. (1997) for the scale-free n = —1 case, since, although the effective n in the CDM 
model does increase from galactic to cluster mass scales, it does not exceed —1 on the cluster scale, so the Tormen et al. 
(1997) results, too, should show a greater degree of central cuspiness than does NFW and they do not. Nor does this explain 
why the simulations of Fukushige and Makino (1998) of a galactic halo, rather than of a cluster halo, also show a density cusp 
steeper than the NFW profile. A second possibility is that, despite the improvement in mass resolution afforded by these new 
simulations, numerical effects are still responsible for a failure to resolve the inner regions of the haloes properly. Splinter et 
al. (1998) for example, have raised a significant question about whether artificial collisional relaxation effects in N-body codes 
used to model collisionless matter have been underestimated, suggesting that previous claims of length resolution below the 
mean interparticle separation are not correct. Further work will be necessary to resolve this matter. 

Regarding the question of the possible difference between the gas and collisionless dark matter distributions, we have 
assumed this difference to be relatively small in our TIS derivation, so that our TIS solution can be applied to either 
component, in proportion to its relative share of the total mean matter density in the unperturbed background universe. This 
was motivated, in part, by the self-similar spherical infall solution of Bertschinger( fl985| ) which supports this assumption if the 
gas has a ratio of specific heats 7 = 5/3. It is also supported by coupled N-body and gas dynamical simulations of the CDM 
model (i.e. those which neglect both radiative cooling and external input of energy to the gas, such as by photoionization and 
supernova explosions). For example, cluster simulations by Eke, Navarro, and Frenk (1998) find that the cumulative gas mass 
fraction measured from the cluster center (in units of the universal mean value over all space) is unity outside the virial radius, 
roughly 0.9 between 0.5r v ir and r v i r and drops below this significantly only in the core. Interestingly, while the density profile 
of the cluster dark matter in those simulations is well-fit by the NFW profile, the gas density profile does depart somewhat in 
the center, where it is found to possess a finite core, and it is better-fit overall by a (3— model like that traditionally used to 
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model cluster X-ray gas as a hydrostatic isothermal sphere, with (3 between 0.7 and 0.75 for the average profiles at different 
epochs, closer to our TIS solution than to the NFW profile of the dark matter. The ratio of the virial radius to the core 
radius of the average cluster density profile at z = 0, in fact, was found to be 20, as compared to our TIS solution in which 
the ratio of truncation radius to core radius is 29.4. (Note: our truncation radius is, moreover, a little higher than the SUS 
approximation for r v ir, so if we replace rt by r v ir = 0.5r m , our ratio drops to 26, even closer). Bryan and Norman(199£) 


report that their simulations of cluster formation in the CDM model by a different numerical technique also support the 
assumption that the cluster gas density profile follows the dark matter density profile throughout most of the cluster volume, 
with differences important only at r/r v n < 0.04. In short, our approximation that the same density profile can be used to 
describe both gas and dark matter in the postcollapse virial equilibrium of a top-hat density perturbation seems justified, 
in view of these simulation results and the uncertainty described above in determining the true nature of the dark matter 
profiles at small radii (i.e. radii as small or smaller than the core radius of our TIS solution). 

We note that, while the dark matter in these simulated haloes in the CDM model is generally found to exhibit velocity 
dispersion which is somewhat radially-biased outside the core, the bias is relatively modest except at radii close to r v n- In 
particular, if f3 an = 1 — v 2 /2u 2 , where vt(v r ) correspond to tangential (radial) velocities, then Eke et al. (1998) report (3 a n < 0.2 
out to r = r v i r /2, rising only to 0.4 at r « r v n. The gas, of course, cannot have an anisotropic temperature but could share 
the radial bias in the sense of there being a net infall or outflow, if it were not hydrostatic. The radial velocity bias of the 
dark matter in the outskirts of the cluster need not be inconsistent with our assumption that the gas follows the dark matter 
as long as it reflects the transition from infall at larger radii to hydrostatic equilibrium at smaller radii for both the gas 
and the dark matter, as already described above to explain the steeper r -3 dependence of the NFW profile at these radii, 
compared to r ~ 2 for the TIS solution. The ratio of dark matter velocity dispersion to gas temperature is reported to be close 
to unity throughout the simulated clusters, except at the very center, where the gas remains isothermal but the dark matter 
“temperature” drops, inside the gas density core radius, in the inner region where the cuspy density of dark matter is claimed 
to occur. In short, there are many indications from these simulations that our TIS solution is a plausible approximation for 
the behaviour of the gas, and also some indication from the similarity between the gas and dark matter outside the gas density 
core that our assumption that gas and dark matter follow the same solution is valid. 

Finally, we return to the question of whether the NFW profile for the dark matter in haloes that form in the CDM model 
is in agreement with the observed properties of dark matter haloes in nature. A full treatment of this question is outside the 
scope of our discussion. We focus, instead on a few recent attempts to answer this question on two widely different scales. 
The first involves dwarf and low-surface-brightness (LSB) late-type galaxies. Such systems are believed to be dark-matter 
dominated, and observation s of t heir rotatio n curv es have been used to probe the density structure of dark matter haloes 
directly. Flores and Primack( [l994|) and Moore( 199-1 ) used rotation curve measurements for several such galax ies and concluded 
that the centra l den sity distribution could not be as cuspy as the r _1 shape of a NFW profile. Burkert( 1995 ) showed that the 
data of Moore( [l994| ) was well fit by a profile of the form 
Po 


p(r) = 


(i + —) 

1 + (—) 2 

V rs ) 

V rs ) 


(107) 


which differs significantly from the NFW profile in that it exhibits a finite core. Kravtsov et al.(1998) have examined this 
question with a somewhat larger observed sample of dwarf and LSB galaxy rotation curves and concluded that a good fit 


to these halo density profiles of the form in equation (106) involves a shallow central profile p oc r r , with 7 = 0.2 — 0.4, 
shallower than the NFW profile. Burkert and Silk( 1997) reexamined the rotation curve of one dwarf galaxy, DD0154, to 
suggest that the discrepancy with the NFW profile could be resolved if a spheroid of dark baryons were added to the CDM 
halo of comparable mass. In short, on galaxy scales, the well-studied systems which are considered to be the most direct 
measure of the dark matter halo are found to be inconsistent with the NFW in requiring either a finite core or a shallow inner 
profile. This is the same sense as the discrepancy between the TIS and NFW profiles. 

On the larger scale of clusters of galaxies, the observations of the X-ray emission by the gas and of gravitational lensing 
by the dark matter provide two independent measures of the distribution of the dark matter in the cluster. Makino, Sasaki 


and Suto(1998) found that the gas density profile of an X-ray cluster, if modeled as an isothermal sphere in hydrostatic 


equilibrium with a dark matter potential well which follows the NFW density distribution, yields an isothermal (3 —profile 
with a finite core, as required to match the observations of X-ray cluster brightness profiles, but with too small a core radius 


relative to the observed ones. Eke et al. (1998) performed SPH simulations of the formation of such clusters in a CDM model 


and claim that their cluster dark matter follows the NFW profile while the gas density profiles are well fit by the isothermal 
(3— model, with a typical core radius for the latter of about 100 h -1 kpc at z = 0, in a good agreement with observations. The 
origin of the difference between these two conclusions is not clear. 

An independent comparison of the dark matter density profiles in observed clusters of galaxies and the NFW profile is 
provided by a recent observation of strong gravitati onal lensing of a background galaxy by the galaxy cluster CL 0024+1654 
at a = 0.39 by Tyson, Kochanski, and Dell’Antonio) 1998 ). These authors construct a high-resolution map of the central mass 
distribution of the cluster and find it to be inconsistent with the singular density profile of NFW. The observed cluster requires 
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a smooth dark matter distribution with a 35/i _1 kpc core and a slope that is slightly shallower than that of an isothermal 
sphere. This is consistent with the TIS solution presented here. 


8.4 Comparison of the predicted virial temperature and mass — radius — temperature scaling laws of the 
minimum-energy TIS solution with numerical simulation results and observations of X-ray clusters 

We have noted above that our minimum-energy TIS solution temperature is larger by a factor of br = 2.16 than the 
postcollapse virial temperature calculated in the SUS approximation, while it is smaller by a factor of 3/2.16 = 1.39 than 
the virial temperature in the limit of a singular isothermal sphere. How does the TIS temperature compare with the results 
of numerical simulations of cosmological structure formation? The clearest comparison available to date is with the results of 
numerical gas dynamics simulations of X-ray cluster formation in the CDM model. It is customary to describe these results by 
comparing the numerical gas temperatures with the so-called “virial temperature” T v n of the clusters. Unfortunately, different 
authors define their T v n differently, often without explicit justification, which makes a direct comparison more difficult. We 
shall attempt to compare the results quoted in the literature uniformly by referring to the definition of virial temperature 
used by Kitayama and Suto(1997), given by 

G[im v M 


H l vir — qV 


3r v 


(108) 


in which r v i r is the virial radius, which different authors have chosen to define in different ways when analyzing their simulations 
(e.g. sometimes r v n is set equal to r 2 oo, the radius within which the average density is 200 times the mean background value), 
and 7 vir is a “fudge factor” which typically ranges between 7 v i r « 1 and 7 v ir « 1.5 and is chosen either to match simulation 
results or according to some analytical model like the SUS approximation. Kitayama and Suto( pT97| ), for example, adopted 
the value 7 vir = 1.2 based upon the gas dynamical simulations of White et al. (1993), while pointing out that other authors 
adopted other values in their analyses, including 7 vir = 1 (Kitayama & Suto 1996), 7 vir = 1.1 ( [Viana and Liddle 1996 ), and 
7 vir = 1-5 (Eke et al. 1996). A complete review of which authors adopt which value of 7 v i r in defining T v \ r is outside the scope 
of this paper. We shall focus, instead, on some recent papers for which simulation results can be used to determine the value 
of 7vir ■ 

We note that the same expression used in equation (|lOS|) can be used to define the virial velocity dispersion, a 2 , of 


the collisionless dark matter component, according to equation (till), if we replace fcsT v i r /(/xm p ) in equation (108) by a 


Henceforth we shall distinguish the 7 vir —values for the gas and the dark matter by the subscripts “T” for the gas and “<r” for 
the dark matter, and drop the subscript “vir” for simplicity. Our assumption in this paper is that the gas and dark matter 
components each follow the same TIS profile with component densities which are simply proportional to the total density at 
each point, and that the two components have the same virial temperature or equivalent virial velocity dispersion a 2 . The 
latter property is sometimes quantified by defining a parameter (3 a -T = (i.e. = 'Ya/'yr )■ We have assumed 

here that (3a-T = 1 is a good approximation. 

A comparison of our TIS solution result, the SUS and SIS approximations, and the results of some of the more recent 
numerical simulation studies of cluster formation in the CDM model is presented in Table For the TIS solution evaluated at 
the boundary £ = </t, = 7 ^ = (6/5 )(&t??) = 1-43 and (3a-r = 1, while the SUS approximation corresponds to 7 t = 7 a = 0.6 

and (3a— t = 1, and the SIS approximation yields = 7 ^ = 1.5 and /3a -t = 1- The numerical gas dynamical simulations 
all report (3a-r ~ 1, with 7 t and 7 ^ ranging between 1.2 and 1.5. The purely collisionless IV-body simulations (without gas 
dynamics) also find 7 between 1.2 and 1.5. These results are in good agreement with the prediction of our TIS solution, 
in rough agreement with the SIS approximation and in clear disagreement with the SUS approximation. At a more detailed 
level, there are some disagreements between the TIS solution and the simulation results and amongst the simulation results 
themselves. For example, Eke et al. (1998) report that the region of hydrostatic and virial equilibrium is most reliably located 
at r ~ r v i r /3 and Evrard et.al. (1996) (hereafter EMN96) note a similar trend, although departures from equilibrium are 
never very great even at r v i r . Similarly, while the clusters are found in all cases to be approximately isothermal, there is some 
tendency for the temperature to fall by less than a factor of two roughly between r v ir/2 and r v i r - Bryan and Norman(199E) 
find a small tendency for (3a-r to decrease with cluster mass, smaller than unity for smaller clusters, but approaching unity 
for the larger clusters, which they attribute to the effect of numerical resolution in better resolving the more massive clusters. 
Finally, Eke et al. (1998) find a small drop in the dark matter cr 2 relative to the gas temperature T near the center, consistent 
with their finding that the central density of the dark matter is more centrally peaked than is that of the gas. Some of 
the differences amongst the simulation results may reflect a weak dependence of those results on the different background 
cosmology models adopted. Nevertheless, within the uncertainties of the simulation results, there is good global agreement 
between the TIS solution presented here and the numerical simulations of cluster formation in the CDM model in terms of 
the temperature - mass - radius virial relation. 

EMN96 have used a large set of simulated clusters to infer a mass-temperature scaling law and a related radius-temperature 
scaling law for the clusters. While the functional dependences of mass and radius on temperature in these scaling laws are 
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Table 2. Comparison of the minimum-energy TIS solution and X-ray cluster simulation results for the 
mass-temperature virial relation 


Source 7 t 7 CT /3 ct _t 


Simulation Results 

Eke, Navarro, and Frenk (1998). 1.5 1.5 1 

Bryan and Norman (1998). 1.2 1.2 1 

Evrard, Metzler, and Navarro (1996)... 1.4 1.4 1 

Crone and Geller (1995). — 1.2 

Cole and Lacy (1996). — ~ 1.5 


Analytical Results 

Uniform Sphere (“SUS”). 0.6 0.6 1 

Singular Isothermal Sphere (“SIS”). 1.5 1.5 1 

Truncated Isothermal Sphere (“TIS”, this paper) ^j. .. 1.43 1.43 1 


a For the TIS solution, 7 v ; r is evaluated at the outer boundary, C = Ct- The values for other radii are plotted 
in Figured, showing that 7 v ; r varies roughly between 1.2 and 1.5 for radii between 0 and rt- 


consistent with the crude expectations of virial equilibrium as measured at a fixed value of the average density contrast inferior 
to some radius r within the cluster relative to the mean background density pb , p/pb, namely, T ~ GM/r ~ ( p/pb)pbr 2 , the 
coefficients of the scaling laws are determined empirically from the numerical gas dynamical simulations of X-ray cluster 
formation in the CDM model. By making a conservative choice, p/pb = 500, of an average density contrast which is generally 
expected to occur at a radius which is significantly smaller than that within which the cluster is well-approximated by 
hydrostatic and virial equilibrium, EMN96 are able to calibrate these scaling laws numerically, with relatively small scatter, 
with the intention of comparing them with X-ray observations of clusters to determine reliable mass estimates for the clusters. 
For our purposes, these numerically - calibrated scaling laws provide an excellent check on our analytical TIS model. If our 
TIS model solution proves to be in good agreement with these scaling laws, then we will not only have confirmed our TIS 
analytical predictions but also have found in our TIS solution a derivation of these numerically - calibrated scaling laws. 

Let rsoo be the radius of the sphere within which the mean density contrast is p/pb = 500, and let A /500 be the mass 
enclosed by this sphere. EMN96 find 


A /500 (T) 
r 5 oo(T) 


/ T \ 3 / 2 

(1.11 ± 0.16) x 10“^) h~'M 0 , 

(l- 24 ± 0.09) 1/3 ft” 1 “PC, 


(109) 

( 110 ) 


where T is an appropriate average temperature for the cluster gas. The simulation results which led to equations (10£) and 
(11C) were analyzed at 2 = 0.04, but scaled to 2 = 0 assuming that they both scale as (1 + z)~ 3 ^ 2 , so there is no information 
explicitly contained in these equations regarding the scaling with redshift. The corresponding scaling laws for a mean density 
contrast p/pb = 200, were presented in Arnaud and Evrard (1998) (with no error bars specified, but the error bars should be 
larger than those at smaller radii; Evrard, private communication) for the same numerical simulation results as in EMN96, 
according to: 


/ r V \ 3/2 

M 200 (T) = 1-45 x 10 15 
r 2 oo(T) = 1-85 (j^) 7 h-'Mpc. 


(111) 

( 112 ) 


We can derive the corresponding relations according to our minimum-energy TIS solution. Let X be the average density 
contrast inside a sphere of radius rx which encompasses a total mass Mx■ Then the minimum-energy TIS solution derived 
here can be used to solve for Mx and rx for the postcollapse virialized object of virial temperature T which results from a 
top-hat which collapses at any z c oil, for any overdensity A', as follows. Let px be the average density inside this sphere of radius 
rx- In terms of pb(t c oil), the mean cosmic background density at the time f C oii for the top-hat perturbation which created the 
postcollapse virialized object described by the minimum-energy TIS solution, the density contrast is X = px / Pb(t co i\) ■ We 
shall continue to refer to the values of these quantities at the outer radius rt by using the subscript “t” as we have in previous 
sections. I 11 terms of the dimensionless radius £ and mass M, the fact that px oc Mx/r\ allows us to write 


MxK 3 x = X 
A/t/Ci 3 Pt / Pb (/coll) ' 


(113) 
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Figure 6. The dimensionless radius » and the corresponding mass Mx , within which the average density with respect to the background 
is X , and the factor 7 v i r > as functions of the density X for our minimum-energy TIS solution. 


Equations (102) and (113) combine to yield an implicit equation for £x and, hence, for Mx, which can be solved 
numerically for any given X , according to 


m = 1871- 2 
M x 


(t )‘ 


The mass Mx can be expressed in terms of the resulting Mx, using equations ([ftj), (^6|), and (JgtJ) , according to 


Mx _ _J_ 

M x ~ 3ttG 


(I) 


-3/2 


M t (1 + ^coll) 


-3/2 H -1 ( kBT 
\nm v 


3/2 


(114) 


(115) 


Similarly, the radius rx can be expressed in terms of the dimensionless radius C,x which is the solution of equation (114) and 
the core radius r o, according to 


rx _ M x 
(x ° k B T M x 


1 / b T \ 
_ 3tt V 5 ) 


-3/2 


Mf (1 + Zcoll) 


-3/2 rr-1 
n 0 


kgT 

/iflip 


1/2 


(116) 


where in expressing ro we used equations (E8I) and (|34[). Equations (114), (115), and (11(:) constitute the final scaling laws for 


the mass and radius within our TIS for the interior spherical region of average density contrast X, in terms of the predicted 
virial temperature T, collapse epoch z co u and Hubble constant Ho, if we use the results of our TIS solution to evaluate br 
and Mt. Plots of C,x and Mx as functions of X, obtained numerically from the TIS solution, are shown on Fig. along with 
a plot of yvir(A') = 2>Mx/C,x which enters equation (108) if we wish to evaluate equation (105) using any radius r < rt and 
replace ill by M(< r) and ?Vir by r. 

above, 


In order to compare directly with the empirical scaling laws of EMN96 in equations (10£), (110), (111), and ( |l 12 
we set X = 500 or X = 200, accordingly, and use our minimum-energy TIS solution summarized in Table [l] to evaluate 
equations (114) - (116). We find £500 = 15.96, M 500 = 37.74, £200 = 24.20, and M 200 = 52.65. The analytical TIS scaling laws 


which result are the following 


/ m \ 3/2 ,/ 11 \ - 3/2 

= ,.llxl0“(—) (1+ (jy 

/ r T 1 \ 1/2 / / // \ — 1/2 

’■“W = 1,24 (lOkev) + 

and 

, rp . 3/2 , \ - 3/2 

M 2 oo(T) = 1-55X10 15 (—) (l + ,con)-^(^) h~'M*, 

/ T \ 1/2 ^ ». x —1/2 

= L 88 (i^) ( 1 +z -^ 3 / 2 (( 4 ) /i_lMpc - 


(117) 

(118) 

(119) 

( 120 ) 
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Table 3. Comparison of the minimum-energy TIS solution and X-ray cluster simulation results for the 
radius-temperature relation. 

X = p/p b Cx M x [rio(X)] BMAr96 0 [rio(X)] T/s {TIS - EMN95)/EMN9e 


200. 24.20 52.65 3.7 3.759 0.016 

250. 21.88 48.59 3.37 3.398 0.008 

500. 15.96 37.74 2.48 2.479 0.0006 

1000.. . 11.56 28.69 1.79 1.796 0.003 

2500.. . 7.32 18.21 1.11 1.137 0.024 


“ Simulation results are from EMN96, Table 5, p. 505, with the value for X = 200 as quoted in Arnaud 
and Evrard (1998). 


Adopting the fully ionized gas value p = 0.59 and the value of 2 co ii = 0 appropriate for a direct comparison with the 
simulation results of EMN96, equations (117), (Ilf), (119), and (120) above become 


M SO o(T) 

r&oo(T) 

and 

A/200 (r) 
T200(T) 


= 1.11 x 10 1 


(loiSv) 

1/2 


3/2 


= 1.24 


(l^v) 


h -1 Mpc. 


= 1.55 x 10 1 


(l5Bv) 

1/2 


3/2 


- 188 (low) 




h 1 Mpc. 


( 121 ) 

( 122 ) 

(123) 

(124) 


These are in remarkably close agreement with the empirical relations of EMN96 in equations (109), (11C), (111), and (112), 
differing from them in radius by less than 2% and in mass by less than 7%, well within the error bars of those empirical 
relations, an excellent check on our TIS solution. In that respect, our TIS solution, based upon the simplifying assumptions of 
spherical symmetry, isothermality and hydrostatic equilibrium, provides an approximate analytical derivation of the numeri¬ 
cally calibrated scaling laws! This close agreement between the predictions of our analytical TIS solution and the numerical 
simulation results which arise as the dynamical outcome of more complicated initial conditions suggests that our simplifying 
assumptions are reasonable. 

A further comparison with the numerical simulation results of EMN96 for the cluster gas is possible in which we compare 
the coefficients of the radius - temperature relation at different radii, corresponding to a range of density contrasts X, derived 
by us according to equation (114) above, with the numerical values reported in Table 5 of EMN96. EMN96 fit the simulation 
results for 2 = 0.04 (divided by (1 + z)~ 3 ^ 2 to yield the scaling law at 2 = 0) and h = 0.5 to a relation of the form 

1/2 


rx(r)™r 10 (X)(Ml_) 


(125) 


A comparison of the analytical TIS predictions for rio(A') in this case, with the tabulated values of EMN96 at X = 
200, 250, 500, 1000, and 2500 is shown in Table |^. The TIS and simulation values are once again in remarkably close agreement, 
the same to within 2.5 per cent at all radii, spanning over an order of magnitude of density contrast, including the canonical 
value of X = 200 so widely used when cosmological simulations are analyzed in comparison with observations of X-ray clusters 
for a wide range of purposes. 

Recently, Hjorth, Oukbir, and van Kampen (1998) have provided a direct comparison between observations of X-ray 
clusters and the theoretical mass-temperature virial relation above, which EMN96 report as a good fit to simulation numerical 
results and which we have derived from our TIS solution. Hjorth et al. (1998) calibrate this relation using observations of 
eight clusters in which the temperatures T are determined from recent X-ray data from the ASCA satellite while the cluster 
masses are determined independently by gravitational lensing measurements. They find that the EMN96 mass-temperature 
relation is in good agreement with their observationally-calibrated relation (i.e. within the observational uncertainties of 
roughly few x 10%). In that sense, we can also claim good agreement between our analytical TIS model mass-temperature 
relation and these observations. 

Regarding the numerically-calibrated radius-temperature relation of EMN96 in equation (11C), which we have reproduced 
exactly here by our analytically derived TIS solution, Mohr and Evrard (1997) have shown that observations of nearby X-ray 
clusters support this relation. Once again, the analytical prediction of the radius-temperature relation by the TIS solution 
presented here is, therefore, consistent with the observations of X-ray clusters. 
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In short, for all its simplicity, the TIS solution derived here appears to predict the virial temperature and the integrated 
mass distribution of the X-ray clusters formed in the CDM model, according to detailed, 3D, numerical simulations of the 
latter, remarkably well. It also matches the mass-temperature and radius-temperature virial relations observed for nearby 
X-ray clusters. We are encouraged by this to think that the minimum energy TIS solution presented here will be quite 
useful as a tool in semi-analytical models of galaxy and cluster formation and as a guide for understanding both the detailed 
numerical simulation results and observations. Further study of the dynamical origin and stability of this equilibrium solution 
and of the sensitivity of the relaxation process which leads to this equilibrium state to initial conditions is warranted. In 
addition, it may be useful to generalize the results presented here to other cosmological background models, such as an open, 
matter-dominated model, or a flat model with a nonzero cosmological constant, in order to describe objects that form at late 
enough times in such models that departures from the Einstein - de Sitter model may arise. We will address this issue in a 
future paper. We expect that, for an open, matter dominated model, the structure of the postcollapse equilibrium object will 
be the same as that derived here for the Einstein - de Sitter case, when expressed in dimensionless units and when referred 
to the quantities at the maximum expansion epoch of the parent tophat. In other words, the density profile will still be as 
plotted in Figure 3, where the core and truncation radii are expressed in units of r m and the density is in units of psus- In 
that case, since the parent top-hat in an open, matter-dominated background universe evolves at a different rate relative to 
the background universe as compared to one in an Einstein-de Sitter model, r m and psus will be different. 
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APPENDIX A: AN ACCURATE ANALYTICAL FITTING FORMULA FOR THE DENSITY PROFILE 
OF A NONSINGULAR ISOTHERMAL SPHERE 


The isothermal sphere solution can be approximated quite accurately by the simple analytical function 
A B 

P ~ a 2 + C 2 ^ b 2 + C 2 ’ 


(Al) 


as suggested by NL97. In what follows, we adopt an approach similar to that of NL97 to solve for the best values of the 
constants in equation (|Al|), different from those in NL97, so as to provide a better fit over the full range of radii required 
to describe the TIS solution derived in this paper. We note that, while the final result of NL97 in the body of their paper 
is correct, several formulae (specifically equations (A9) and (All)) in their Appendix A2 are incorrect and thus the method 
should be exercised with care. 

The solution of equation (^9j) has a Taylor expansion at C = 0: 


- _ i 1 a 2 . 1 a 4 61 6 

P ~ 1_ 6 C + 45 C ~ 22680 C + 


629 


2041200 


C 8 + 0(C 10 ). 


(A2) 


Depending on our requirements on the solution, we can find diffe rent sets of parameters for the best fit. If, 
require that the first three terms of the Taylor series of equation (© at C = 0 are the same as in equation 
has the right behaviour at infinity, 


for e xample, we 
(A2) and that p 


2/C 2 as C ■ 


then we obtain the simple approximation of equation (A2) with the following parameters 
(A, a 2 , B, b 2 ) NLA = (50,10, 48,12). 


(AS) 

(A4) 


This approximation is fairly accurate for small radii (fractional error less that 2 per cent for C < 5) and for very large radii 
(error < 2 per cent for C > 100). However, in the intermediate interval, which we are most interested in, the fractional error 
of this approximation can be as large as 24 per cent. To remedy this problem, NL97 offers the following choice of parameters: 

(A,a 2 ,B,b 2 ) NL ,2 = (44.34,10.16,42.61,12.66). (A5) 
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Figure Al. Fractional error (p exa ct - Papprox)/Pexact 


of the approximation of equations (|Al|) and ( 


This latter approximation is an extremely good fit up to £ = 10, but deteriorates after that, having fractional errors of 5-8 per 
cent for £ = 15 — 40. Hence, we have used their method, instead, to find an approximation which better suits our purposes, 


by requiring the first three terms of the Taylor series to match those in equation (A2), and by imposing an exact match at 


Q = 29.4, which produces the result quoted in our equation (103) in the main text, namely 
(A, a 2 ,B, b 2 ) T is = (21.38, 9.08,19.81,14.62). 


(A6) 


While we cannot guarantee that this is the best possible approximation of this type, it is accurate to within 3 per cent for 


£ < 40 (Fig Al), while for most of the interval the error is much smaller, which should be sufficient for many purposes where 


using a simple functional form is more important than obtaining a high precision. 
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